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ABSTRACT 

Atmospheric temperature structure fluctuations quantified by 

degrade the spatial and temporal coherence of electromagnetic and 
acoustic waves propagating in the atmosphere. A computer controlled 
atmospheric echosounder developed at the Naval Postgraduate School 
measures a time averaged Cy^ profile of the lower atmosphere. Assigning 

the proper Cy^ values to the backscattered return signals depends on an 

accurate calibration of this instrument. Calibration involves 
determination of the product of the echosounder's transmission efficiency 
and reception efficiency E^. 

This thesis provides a preliminary investigation of a calibration 
process using pulsed acoustic energy backscattered from hard spheres. 
Supporting software calculates the desired product based on an 

assumption of echosounder efficiency reciprocity. Results of the 
calibration process investigation indicate this assumption may be invalid. 
The results also indicate the software performs as intended and that the 
proposed calibration method possesses sufficient merit to warrant further 
development. 
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I. 1NT.BQDU.CI1.QN 



The variational effect of wind shears, convection, and temperature 
gradients on atmospheric density distributions characterize atmospheric 
turbulence. The temperature structure parameter, Cj'^, and the index of 

refraction structure parameter, quantify two density distributions 

that arise in acoustic and electromagnetic propagation. [Ref. 1] 

The density distributions quantified by and induce phase 

fluctuations in propagating electromagnetic waves. The density 
distributions of atmospheric turbulence degrade the spatial and temporal 
coherences of initially coherent electromagnetic radiation. Similar, more 
pronounced effects occur in the propagation of acoustic waves. [Ref. 1] 

A number of techniques for correcting these problems have been 
proposed, particularly with respect to electromagnetic propagation. These 
corrective techniques include the use of modular mirrors [Ref. 2], optical 
phase conjugation [Refs. 2 and 3], optical and digital signal processing 
[Ref. 4], and simple avoidance of turbulent air masses. The appropriate 
timing of transmissions or suitable location of transmission and reception 



facilities accomplishes avoidance (Refs. 5 and 6], These techniques can 
be applied alone or in combinations. 

A knowledge of the history of the density or refractive structure 
along a proposed transmission path is required to effectively employ a 
number of these compensatory techniques. Time-averaged or 

profiles charactize this history [Refs. 1 and 7], 

Wroblewski and Weingartner developed and Moxcey refined a high 
resolution, computer controlled acoustic sounder, or echosounder, system. 
Their computer controlled echosounder system uses a transmitted pulse of 
acoustic energy as a probe of atmospheric structure. The interaction of 
the acoustic wave packet with atmospheric turbulence is used to develop a 
time-averaged profile for a short spatial range [Refs. 5, 6 and 8]. The 

acoustic sounder may also be used to verify certain long range optical 
measurements of turbulence [Ref. 6]. 

Development of the profile uses a expression developed from 

the echosonde equation summarized by Neff [Ref 7] and the empirical 
acoustic backscatter cross section per unit volume expression of Tatarski 
[Ref. 9] [Refs. 5 and 6]. The accuracy of the absolute values of this C.^.^ 

profile depends on the calibration of the acoustic antenna parameters 



2 



contained in the expression. The antenna parameters of interest are 
the efficiency of conversion of electrical to acoustic power, and E^, the 

efficiency of conversion of acoustic power to electrical power. 

In this thesis acoustic energy backscattered off acoustically hard 
spheres was used to determine E^, and E^. Acoustic range equations were 

software implemented to perform the necessary calculations. Data and 
calculational results provided first order verification of the calibration 
process. 
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II THEORETICAL BACKGROUND 



A. REVIEW 

Wroblewski and Weingartner applied the echosonde equation, 

P.= E,[p,E.][e-^“^][^o„(R,0][A 6]_ (1) 

summarized by Neff [Ref. 7] and the acoustic backscattering cross 

section per unit volume expression given by Tatarski [Ref. 9] as 

1 / 

(Jp(R,f) = 0.0039 k 

to develop the volume averaged expression 

C " = ^ 2 _ _L 

^ 0.0039 E^E^ CT AG 



(2) 



Pt 



(3) 



where 

• is the electrical power of the reflected signal returned to the 
acoustical antenna array, 

• is the efficiency of conversion of the returned acoustic power to 

• is the electrical power supplied to the acoustical array, 

• is the efficiency of conversion of P^ to acoustic power, 

• a is the average attenuation per unit distance along the transmission 
path, 

• R is the range from the array to the target or backscattering volume. 
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• c is the average speed of sound along the transmission path, 

• T is the length in time of the transmitted acoustic pulse, 

• A is the antenna’s aperture area, 

• G is the antenna's effective aperture factor, 

• k is the wavenumber of the incident acoustic energy, and 

• Tq is the average temperature (in degrees Kelvin) along the 

transmission path. 

[Refs. 5 and 6] 

The conditions at the time and location of each application determine 
the value of the terms T^, k, c. x, P^, P^, R and a in the expression. 

The terms A and G depend on the antenna design. 

Adverse field effects such as dust, debris and vibration loosened 
connections cause the remaining two terms, and E^, to change over the 

life of the system. This chapter develops the theoretical expressions 
supporting a determination of E^ and E^^. 



B. ACOUSTIC RANGE EQUATIONS AND EFFICIENCIES 

In the following two sections the acoustic range equation of Neff 
[Ref. 7] is re-developed for targets of small area relative to the 
cross-section of the ensonifying field. Neff [Ref. 7], Skolnik [Ref. 10], and 
Probert -Jones [Ref. 1 1] were used for guidance. 
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1. One Way.Range Equation and EtGo 



As defined previously, the electrical power supplied to the 
acoustic array is P^. It is converted at an efficiency to an emitted 

acoustic power, 

Pa = Pt E, . (4) 

Pg is spread over a 4tt solid angle. 

A far field range R is defined by the condition 

X/2 » ILI - IRl (5) 

for the geometry in Figure 1. A is the wavelength of the emitted acoustic 
wave. 




Array 



Figure 1; Calibration (Far Field) Geometry 
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The acoustic intensity, 



I 3 = Pg / 4 tir2, ( 6 ) 

describes the emitted acoustic power density of an isotropic emitter at a 
far field range R. Ig is in units of [W/m^]. 

If the emitted power is directional, the intensity Ig is modulated 
by a geometrical gain, 

G(O) = GglF(O)|2 = Gg|F(0,8)|2 (7) 

Note this is not the same 6 found in equation ( I ). Now, 

|F(0,0)|2=1, (8) 

and G^ is thefore the centerline gain of the acoustic array. [Ref. 1 1] The 

centerline direction (0,0) is taken to be normal to the face of the acoustic 
array. The acoustic intensity Ig is now written 

U = — r GiO) (9) 

In a real propagating medium intensity attenuation occurs along 
the transmission path. The radial term e*^^ accounts for this 
attenuation when the atmospheric path is assumed to be reasonably 
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homogeneous. The acoustic intensity is then 



_ G(0) 



4ttR^ 



( 10 ) 



A target of cross sectional area intercepts a power 



*^tgl ’a ^Igt 



( 11 ) 



Expanding this expression with equation ( 1 0) yields 






( 12 ) 



For a target of cross sectional area A^^^ with a sufficiently 
small solid angle Aj_g^ /R^ little change in intensity occurs about the 

direction 0. Equation ( 1 2) is now arranged as a one way range 
equation along the centerline direction (0,0); 

p - r^t^t e‘“''G .^1^ 03) 

- L 4« " p. , • 

igt 



Isolating E^G^ in the above expression yields 



= 



4xf R 



P 

OCR tgt 
'^tgt 



(14) 
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Electrical power is expressed as 



Vrms ^ 

Z Z 



( 15 ) 



where 



• is the empirical root mean square voltage, 

• V is V and 

• Z is the electrical impedance. 



E^Gq is now written 



EtGo 



= 4rf R 



(V ms)^ 



P 

OCR tqt 



‘tgt 



(16) 



The ratio is defined as the acoustic intensity 1^ 

received by the target. Therefore, 



= 4rfR 



OCR 



ms)^ 



(17) 



2. Two Way Range Equation and Er 

For an acoustically hard target the total power reflected is 

The power reflected from the target surface has a normalized, 
directional intensity distribution described by the differential scattering 
cross section, cr. 



9 



The direction from the target to the array is the backscattered 



direction. The differential scattering cross section in the backscattered 
direction is called the normalized backscattered cross section, is 

defined as the power reflected toward the source per unit solid angle, 
normalized by the incident power density over a 4 tt solid angle. 

The backscattered intensity from the target is 



An array of aperture area A at range R intercepts an attenuated, received 
acoustic power. 



Substituting equation ( 1 8) into equation ( 1 9) yields 

P = (T, e-“" A (20) 

4ffR^ 

is converted to returned electrical power, P^, at a return 

la r 

efficiency, E^; 



P, 



ra tgt 



A e . 



(19) 




( 21 ) 
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Since electrical power is expressed as 



P 



Vrms ^ 
Z 




(15) 



electrical power representing the returned acoustic signal from a 
target placed along an emitter's direction (0,0) is written as 

(Vms)r _ (Vms)t ^ 1 (22) 

Zr “ Zt L 4TT r2 Jl 4tT r 2 • 



For a passive, linear acoustic array circuit the concept of 
impedance reciprocity is invoked [Ref. 12] to yield a two way range 
equation for scattering from a small target; 

1 







2 ^ 
(4xf) 




(23) 



This is written to isolate E^.; 

. _ (V»^)r (4-nR^f j_ . 

'■ ■ (V^), E,G„ A A,,, 0^ - 

A part of the return signal's power is unwanted noise. 
Subtracting noise power from the returned acoustic power yields 



Er 



(Vms)r-(Vm5)n (4T)R^)2 

E.G. 



(Vm.)t 



t'-'O 



2<XR 

_e 

A A 




(25) 



As return voltages are quite small they are analyzed after 



electrical amplification. Calling the electrical amplification a gain 



G„ allows introduction of the final modification to E : 



c - ( ms)r~ ( ms)n (4T) ^2 <X R 

(Vms)i ^ 



^ (26) 



C. ATTENUATION AND RANGE DETERMINATION 



Employing equations (17) and (26) in the calibration process requires 
calculating the attenuation a and the range R from the array to the 
target. 

1. Attenuation 



The calculated attenuation cc for an atmospheric propagation 
path is assumed to be the combined result of molecular and classical 
absorption in the atmosphere. The classical absorption oc^, is attributed 

to heat conduction and viscous effects while the larger molecular 
absorption is attributed to the excitation of internal energy modes of 

atmospheric gases by the propagating sound energy. [Ref. 7] 

Information from Neff [Ref. 7], Businger [Ref. 13], Neiburger 
[Ref. 14], and Fuller [Ref. 15] is applied to make a reasonable determination 
of cc. The determination is made using the average measured temperature 
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in degrees Celsius along the propagation path, the measured 

atmospheric pressure P in millibars, and the measured relative humidity 
Rh in percent. 

The molecular attenuation coefficient, is expressed by the 
empirical relationship 



*^mol 



*^rri8x 

304.8 







1 + f 



ratio 



vn 



( 27 ) 



“mol units of [-dB/m]. 

Now, 

= 0,0078 f^ e7.77(i-i/T*) (28) 

is the maximum absorption at 

(to* 6600 ( 1 OO-p- ) ♦ 44,400 ( 1 00 f-f ) P* (jg) 
= (T*)°® . 

f is in units of [Hz]. 

m 



Also, 



f’ . Eg Rh 
P P - Egd + Rh) 



( 30 ) 



is the mole ratio of water vapor in the atmosphere. The water vapor 
pressure e is expressed in millibars. 
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(31) 



E3= 10^^-^-2353/TJ 

expresses the saturation vapor pressure E^. is the average absolute 

temperature along the transmission path, as defined in Section A. 

T* = (1.8T^+ 492)/ 519, (32) 

and 



P^ = P/ 1014. 



(33) 



Finally, 



Vatio ^ 



(34) 



where f is the operating frequency of the echosounder in [Hz], [Refs. 7, 
13, 14 and 15] 

The classical attenuation coefficient, cc ., is approximated by 



a^j = 1.74^ 10‘^° (f)2 . 



(35) 



a^., is in units of I-dB/m], [Ref. 7] 



For a total attenuation 



«T = «cl " «mol 



(36) 



and the relationship 
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- a-j-R = 10 log e 



-aR 



(37) 



the average attenuation per meter is finally written as 



oc , + oc , 

OC = in 10 



(38) 



Equation (37) is derived from unit definitions and equation one of Neff 



[Ref. 7], 



2. R angfi 

The coordinate system for the backscattered cross section 
determination uses the center of the target spheres as the coordinate 
system origin. Therefore, the radius a of a target sphere and the time of 
return t^ of an acoustic signal from the acoustic echosounder are used to 

determine the range R from the target to the echosounder; 

R = . (39) 



c is the average speed of sound along the transmission path, and 

V “ deceive ” ^transmit • 

For the range of temperatures expected in normal echosounder operation 
the speed of sound in meter per seconds in dry air is 

- 20.05 /f; . (41) 
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For moist air c becomes 



^ ■ Vy<' (42) 

where e /9 is defined in equation (30). [Ref. 1] 

D. BACKSCATTERED CROSS SECTION DETERMINATION 

Employing equations (17) and (26) in the calibration process also 
requires calculating g-j^, the backscattered cross section. Recall from 

Paragraph B.2 that the backscattered cross section is the differential 
scattering cross section in the direction of the energy source. 

1. Initial Conditions and Assumptions 

In determining the differential scattering cross section <r it is 
assumed the incident wave state is unchanged after scattering. That is, 
for a stationary, acoustically hard target it is assumed 

^ ” ^incident ~ ^scattered ’ 

and that incident k has good definition, i.e. the incident wave can be 
considered monochromatic [Ref. 16]. 

Selection of acoustically hard spheres as targets further 
simplifies the determination of the differential scattering cross section 
solution. This target selection sets the boundary condition at the target 
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surface as the Neumann boundary condition 



dn 



0 . 



( 44 ) 



1 / is the velocity potential defined by the acoustic propagation velocity c 
as 



c = 7 f-' . (45) 

n lies along the direction normal to the target surface. Also, 

P = - Po if 

for the pressure p and rest density of the propagating medium. [Ref. 16]. 



V^2 
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Figure 2; Incident Wave and Target Description 



The target sphere is aligned along the centroid of the main 
acoustic lobe of the acoustic arrray echosounder, that is, along the 
direction (0,0), normal to the face of the array. The target sphere is 
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selected to be sufficiently small with respect to incident beam 
dimensions so that an essentially uniform intensity is intercepted by 
the sphere [Ref. 17], Since the target intercepts a very small proportion 
of the centroid of the transmitted sounding lobe no accommodation for 
sounding lobe geometry over the target’s angular area is made. 

Acoustic gain is simply taken as Gq. 

Target placement is also at a far field range. This allows 
spherical incident waves to be treated as planar incident waves and 
places the target outside the range associated with the ring time of the 
echosounder [Ref. 8]. Placement of the target is also close enough to 
the array to support an assumption of transmission medium homogeneity 
along the transmission path. 

Finally, the potential describing the target falls off faster than 
1/r. The describing potential for an acoustically hard sphere with 
radius a is 



V(r) = 



r i a 



- 0 r i a 

Thus, for r > a V(r) falls off faster than 1/r. 



( 47 ) 
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Note that with the exception of the specific surface boundary 
condition for acoustic waves the above conditions and assumptions closely 
correspond to those for plane electromagnetic waves incident on a 
perfectly conducting sphere, and for perfectly elastic nuclear scattering 
in a center of mass coordinate system. [Refs. 16, 18, 19 and 20] 

2. The Incident Wave 

The asymptotic form of an incident plane wave approaching a 
scattering target but beyond the influence of the target's potential V(r) is 
the undisturbed, or free, plane wave. In the spherical polar coordinates 
(r, 0 ,e) illustrated in Figure 3 a free, monochromatic plane wave of 
unit amplitude propagating in the direction 0 ^^ given by ( 00 , 8 ^,) is 
f(r,t) = e ^ sin(0o)cos( 8 - Sq)] g-iCDt 

[Ref. 16] 




Figure 3: Spherical Polar Coordinates 
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For a wave of unit amplitude incident along the z axis , 0^ = 0 
and the free wave f(r,t) is presented as 

|(r t) = e il^rcos(0) g-M ^ g ikz g-M 

Note this is a solution with axial symmetry for the free wave equation 
[Ref. 18], 

The wave equation describing the behavior of the undisturbed or 
free monochromatic wave |(r,t) is 

( = 0 . (50) 

at" 



c is the speed of propagation of the free wave. If the wave 
function f(r,t) has a simple harmonic time dependence, i.e., 

f(r,t) = e-*“* 1(c) , (51) 

then for a wave number 

k = (o/c = 2 ti/A (52) 



the free wave equation becomes the Helmholtz equation 

[v2.k2]f(r,t) = 0 



[Ref. 211 



(53) 
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The stationary solution j(r) of the free wave equation can be 
written as the product 

|(r) = ||(q| )| 2 (Q 2 ) |3(Q3) (54) 

for orthogonal curvilinear coordinates q^ d 2 , Q 3 [Ref. 16], Recall an 

undisturbed plane wave has a cylindrical symmetry about the axis of 
propagation [Ref. 22]. Using the spherical polar coordinates 

q,= r, P 2 = 0 , q 3 = 0 , (55) 



the stationary solution describing an undisturbed plane wave propagating 
along the z axis is 

f(r) = f,(r)| 2 ( 0 ) . (56) 

Substituting f(r) into the free wave equation expressed in 
spherical coordinates and separating the coordinate dependencies yields 
the separated free wave equation, 



= - - — ^ — (sin 0 — f ). (57) 

flh I f2Sin0 d0 d0 ^ 

Using a separation constant of 1(1+ 1 ) the equation for the angular 
dependency of the free wave becomes 



7 —! ^ (sin 0 -^ f.) 1(1 + 0 = 0 (55) 

Jo sin 0 d 0 d .0 ^ 
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Letting 



cos(0) = X (59) 

in the angular equation yields 

[(1 - + 1(1 + 1)] = 0 . (60) 

'• dx^ dx ^ 

The previous equation is known as the Legendre differential 
equation [Ref. 21], The most useful solutions to the Legendre differential 
equation are the Legendre polynomials, designated by Pj(x). They can be 

generated by the recursion relationships 

Po(x)=1, (61a) 

P,(x) = x (61b) 



and 



P,(x) = [2 - 1 /I] X P,_ , (x) - [ 1 - 1 /I] P,_ 2 (x) (6 1 c) 



[Ref. 21], Therefore, 



|2 ](0) = P,(cos 0) . 



(62) 



From equation (57), the equation describing the radial 



dependency of the free wave is 

■ dr"^ r dr 



* = ( 63 ) 
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Changing scale to the variable s = kr yields the differential equation 




ZA. 
s ds 



[1 



l(Hl) 

s2 






(64) 



[Ref. 17], 

This differential equation has the spherical Bessel function j,(s) 
with amplitude A, as a solution; 

I, ,(s) = A,j,(s) (65) 

[Ref. 18]. 



The stationary solution of the free wave equation in spherical 
coordinates is now written as the partial wave 

|,(r) = P,(cos 0 ) [A, j,(kr)] . (66) 

The most general form of this free wave equation solution is a Fourier- 
Bessel series, 

OO 

|(r, 0 ) = I^A, j,(kr) P,(cos 0 ) (67) 



[Refs. 17 and 21]. 

A plane, free wave of unit amplitude propagating in the direction 
Oq= 0 is then 

e'^ = 2|.|A, j,(kr) P,(cos 0 ) . (68) 
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The norm of the Legendre polynomials, the asymptotic expression 



j ](kr) ~ 



sin (kr - ^ 0 
kr 



(69) 



and the process detailed in Elton [Ref. 18] are used to determine Aj; 

A, = iM21+n. (70) 

Therefore, 

= ^2^i'[21 +1] j|(kr) P^(cos 0 ) (71) 

[Refs. 18 and 21], 

3. Detection Conditions and Assumptions 

As the dimensions of the incident acoustic wave packet are 
greater than the size of the target it is expected that a portion of the 
incident wave packet continues to propagate past the target. The total or 
disturbed wave description is then a superposition of the incident and 
scattered wave descriptions. [Ref. 17] 

The acoustic array serves as both transmitter and receiver, or 
detector. This keeps the detector outside the path of continued 
propagation of the incident wave packet. [Ref. 17] Therefore, only the 
scattered wave description is of interest in the calibration of the array. 
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The placement of the target at a far field range simultaneously 
places the detector in the region of asymptotic behavior of the scattered 
wave. 



4. The Disturbed Wave 

Once an incident wave described by the velocity potential l/(r) 
enters the region of influence of a scattering target's spherical potential 
V(r) the wave equation becomes the disturbed wave equation 

[ + [k2- v(r)]} ^"(r) = 0 (72) 



[Ref. 1 7], Comparison of the form of the disturbed wave equation to the 
free wave equation gives an expectation of a solution of the form 

l/(r)= l/^(r) 1 / 2 ( 0 ) . (73) 



Substituting //(r) into the disturbed wave equation and 
separating variables yields 



L A (siri0^)K = r^[k^- .(74) 

1/2 sin 0 d 0 d .0 4 '■ l/^ dr dr 1 



The angular equation for the disturbed wave is identical in form 
to the angular equation for the free wave and independent of the Neumann 
boundary condition d I// dr = 0. Therefore, 

1^2 1 ( 0 ) = Pi(cos 0 ) . (75) 
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Using the previous substitution of s = kr yields the now modified 



Bessel equation as the radial equation for the disturbed wave; 

{-^ + [ 1 - r 0 (76) 

^ ds^ s ds s2 ■'-II- 

Following the solution development in Elton [Ref. 18] and 
applying the Neumann boundary condition at r = a leads to a solution 
//^ ^^(r)= iU2W l] [j,(kr) - a^‘ h/’Hkr)] . (77) 



The relationship 



a, = i e"'^i sin 8, = — f-n (78) 

' ' h,''> (ka! 

is defined from Morse [Ref. 17] and Bowman [Ref. 16]. Additionally, 

• the j,(kr) are the spherical Bessel functions of the first kind, 
argument kr. and order 1. 

• hj^i^kr) are the Hankel functions of argument kr and order 1. 

• the 8| are the phase shifts of the usual partial wave description for 
spherical scattering solutions [Refs. 17 and 20]. 

• the j,'(ka) are the first derivatives of the j,(kr). taken with respect 
to the argument kr and evaluated at r = a. 

• h,^i^ (ka) = j,'(ka) i y,'(ka). and (79) 

• the y,'(ka) are the first derivatives of the y|(kr). the spherical Bessel 

functions of the second kind with argument kr and order 1, taken with 
respect to the argument kr and evaluated at r = a. 



26 



The genera] series solution for the total disturbed wave i/{r) is then 

l''{r)= S i’[21+ l] P,(cos 0 ) [j,(kr) - a,' h/ ^Hkr)] . (80) 

1=0 ' ' ' ' 

5. The Scattered Wave and Its Amplitude 

The general asymptotic form for a disturbed wave Kr,0) 

resulting from a plane wave |/j(r,0) incident along Oq = 0 and scattered by 
a spherical potential V(r) is 

|/(r,0) = e’*^^ + s(0) e’^^^/r (81) 

[Ref. 18], This description of the asymptotic disturbed wave is a 
superposition of the incident plane wave e'*^^ and a spherical, scattered 
wave s(0)e''^''/r. The term s(0) is called the scattering amplitude 
[Ref's. 1 8, 1 9 and 22], 

Now, recalling the series equivalence for the plane wave incident 
along 0^ in equation (71), and comparing the general series solution for 

the disturbed wave in equation (80) to its asymptotic form in equation 
(81) gives 

s(0)e’^''/r = - 5 iU2 W 1 ] P,(cos 0 ) a,’ h/ ’ )(kr) . (82) 

1=0 ' ' ' 
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This is a series expression for the scattered wave. Therefore, 



s(0) = ^ - 7- - 1 ? [21 + 1] P,(cos 0) a,' (83) 

k ]=o ' ' 

[Refs. 16 and 18], 

6. The Backscattered Cross Section 

The scattering amplitude s(0) is related to the differential 
scattering cross section a-(O) = by 

<7(0) = |s(0)|2 (84) 

[Refs. 18, 19 and 22], 

Applying the series expression for the scattering amplitude in 
equation (83) to the previous equation yields 

(tA. . =<r(0)= -j IS i [21 + l]a;P,(cos0)|2. (85) 

•■S'- 1^^ 1=0 ' ' 

For the backscattered direction 0 = Tf this gives a backscattered cross 
section of 

4 TT 00 , 1 , n 

[Refs. 16 and 17]. For the spherical target = na^, so 

A 00 - -1.0 

(7. = -^ 1 I [2M i]a, p. (86b) 

^ (ka)^ 1=0 1 
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It should be noted a plane wave of unit amplitude incident from 



the positive z direction and backscattered along 0 = 0 would introduce a 
factor of (-1)' into the previous sum (Ref. 16], 

E. SPHERICAL BESSEL FUNCTIONS AND THEIR DERIVATIVES 

Implementation of the previous equation necessitates calculation of 
the The recursion relations 



jpUa) = [-yka) + cos(ka)]/ka , 


(87a) 


j,‘(ka) = {-[1 + l][j,(ka)]/kal + j,_|(ka) , 


(87b) 


yQ(ka) = [-yo(ka) + sin(ka)]/ka 


(87c) 


y,'(ka) = [-[1 + 1 ][y ,(ka)]/ka} + y ,(ka) 


(87d) 



determine the j,'(ka) and y,' (ka) needed for the calculation of aj' in 
equation (78) (Refs. 21 and 23], 

Determination of the regular and irregular spherical Bessel functions 
at r = a supports the use of the recursion relations in equations (87). As 
ka is real the irregular spherical Bessel functions are stable enough to use 
the recursion relations 
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y^do) = -[cos(ka)]/ka . 



(88a) 



y j(ka) = [y^Cka) - sin(ka)]/ka 



(88b) 



and 



y,(ka) = {[21 - l][y,_,(ka)]/ka) - y ,_2(ka) (88c) 



in their determination [Ref's. 23 and 24 ]. 

Because the size of the argument of the regular spherical Bessel 
functions may be very small the method of negative continued fractions is 
appropriate to the fast and accurate determination of a converging 
solution to these functions [Ref.24j. The regular spherical Bessel function 
reciprocals are therefore determined by the negative continued fraction 
calculation 




Now, 



b^^^ = 2(1 + m - l)/ka 



(90) 
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and 






(91) 




1 



b 



1.2 b 



(Ref. 25], 



While the regular spherical Bessel function reciprocal calculation 
of equation (89) would appear to involve an infinite number of terms, the 
calculation can in fact be terminated when the m^ numerator term equals 
the m^ denominator term to the desired precision for the spherical Bessel 
function's solution. [Ref's. 24 and 25] 
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III. CALIBRATION PROCESS AND 50FTWARF 



Small quantities of data were collected in the anechoic chamber at 
the Naval Postgraduate School. This data aided first order verification of 
the calibration process and its supporting software. A brief description 
of the anechoic chamber is provided in Appendix A. 

The calibration process and data development occured in two stages. 
The first stage used equation (26) and two way propagation path 
measurements to determine (t^^E^Ej^Gq. The second stage determined Ej^Gq 

using equation ( 1 7) and one way propagation path measurements. 

E was then determined from <r. E„E.G„, E.G« and the calculated <r. of 

r brlO'lO d 

equation (86b). 

Because reciprocity was assumed for E^, and E^ 

= E^2 (92) 

provided the desired calibration information. 

Supporting software was written to perform calculations required for 
calibration The next section describes the development of this software. 
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A, CALIBRATION SOFTWARE DEVELOPMENT 



Calibration software was written in BASIC. It was used with a 
BASIC 4.0 compiler on a Hewlett Packard 3 1 0 computer. Computer 
performance was enhanced by an Infotek floating point processor. 

Notation in the programs sometimes differs from that in the thesis 
text as a result of H.P. keyboard limitations. Variable definitions are 
located at the beginning of subprograms or modules. These definitions 
are provided throughout a program as an aid in determining program to text 
variable equivalence. 

The main effort of software development was calculation of 
backscattered cross sections. This was done with equation (86b). 
Information available for summation loop verification was differential 
scattering cross section information. Differential scattering section 
software was developed first. 

1. Differential Scattering Cross Section Software Development 
The differential scattering cross sections a were calculated 
using equation (85) and equation (78). 

To support equation (78) employment the subprograms 
calculating the values of the Legendre polynomials and the values of the 
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regular and irregular spherical Bessel functions were developed first 
Each subprogram was developed and tested separately in test programs. 
Each test program was designed to provide a tabular output for 
verification against values found in Abramowitz and Stegun [Ref. 23], 

The Legendre values were calculated using equations (61 ). The 
irregular spherical Bessel functions were calculated using equations (88). 
The regular spherical Bessel functions were calculated using equations 
(89), (90) and (91). 

Following verification these subprograms were assembled with 
subprograms calculating the derivatives of the regular and irregular 
spherical Bessel functions, the normalized differential scattering cross 
section c, the scattering modulus of Bowman [Ref. 16], and with three 
graphics subprograms. This assembly created the differential scattering 
cross section program found in Appendix B. 

The values of the regular and irregular spherical Bessel function 
derivatives at r=a were calculated using equations (87). 

Differential scattering cross sections were calculated using 
equation (85). The number of orders 1 required for the summation 
calculation in equation (85) was determined by generating tabular outputs 
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and determining the number of orders required for numerical stability to 
occur to 12 decimal places. This stability process was based on the 
concept of the phase shift 6, of each outgoing partial scattered wave 

becoming progressively smaller with increasing 1 until (equation (78)) 

became negligible for the precision desired [Ref's. 18, 19, 20 and 22]. 

Tabular outputs covered a ka value range of one to thirty at 
varying scattering angles, including 0 = Tf. This ka value range selection 
was based on the expected maximum target size to be encountered during 
actual echosounder employment. 

Twenty-one orders greater than the integer value of the 
argument ka were determined to be required for equation (85)'s 
summation. The number of orders required was greater than the usually 
cited relationship 

ka~ Vl(l + 1) (93) 

gives. It should be noted that the precision criteria of (r was greater than 
that normally encountered. (Refs. 18, 19, 20 and 22] 

The graphics subprograms of the differential scattering 
cross section program provided outputs for comparisons to Figure 80 in 
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Morse [Ref. 17] and Figure 10. 10 of Bowman [Ref. 16]. These comparisons 
were used to verify the functioning and output of the differential 
scattering cross section program. 

The graphics subprograms used to generate Figure 4 provide 
linear polar plots of the differential scattering cross section ff(0) 
normalized by the cross sectional area of the target sphere. These 

linear polar plots show the normalized differential scattering cross 
section of acoustically hard spheres ensonified by plane acoustic waves 
incident from the left, i.e. from 0 = tt. 

The first view of each sheet in Figure 4 shows the entire 
normalized differential scattering cross section pattern, including the 
forward lobe, or "shadow zone", along 0 = 0. The second view of each 
sheet is a magnified plot showing the scattering detail at a scale 
comparable to the normalized backscattered cross section value. Note 
the two views are identical on Sheet 1. The need for magnification with 
increasing ka becomes apparent on subsequent sheets. 
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Pl/2 
. t 




C - 343 M/5. 

F » 5000 HZ. 

A » ,0109180290951 M. 

K«A - 1 . 

NORMALIZED 8CKSCTR - .680073464081 



Figure 4, Sheet 1 of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 80 of Morse iRef. 17], 
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pi/a 



3 : 



2 * 




3»P]/2 

P]/2 

t 

i 

1 




3*PI/2 



C “ 343 h/S. 

F - 5000 HZ. 

A = .0327540872883 M. 

K*A - 3 . 

NORMALIZED BCKSCTR - .574927362173 



Figure 4, Sheet 2 of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 80 of Morse [Ref. 17], 
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p ]/2 




i 

6 




I 

4 



3*PI/E 

PI/2 




C “ 343 M/S. 

F - 5000 HZ. 

A «= .0545901454805 M. 

K • A " 5 . 

NORMALIZED BCKSCTR - 1.08120386888 



Figure 4, Sheet 3 of 3; Differential Scattering Cross Section 
Program Output for Comparison to Figure 80 of Morse [Ref. 1 7], 
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Figure 4 was used for comparison to Figure 80 of Morse [Ref. 17] 
The radial scale of Morse was not graduated so verification consisted of 
a check for the correct number and approximate angular location of 
minimums and maximums in the differential scattering cross section 
patterns. 

The graphics subprograms used to geherate Figure 5 provide 
semi-logarithmic plots of the scattering modulus 5 of Bowman [Ref. 16]. 

J = k s( 0 ) / (94) 

is the scattering modulus 5 of Bowman [Ref. 16]. s(o) is the scattering 
amplitude s(0), as defined by equation (83). with a factor of (-1 )' 
introduced into the sum. that is. 

s( 0 )= ^ i z (-1)^21 + 1) P/cos 0 ) a ■ (95) 

k 1=0 ' ' 

s(0) is normally used to calculate a-(O). per equation (84). The scattering 
modulus 5 of Bowman [Ref. 16] could also be used to calculate <r(0). Note 
the introductioh of (- 1 )’ into the summation corresponds to the case of an 
acoustically hard sphere ensonified by a plane acoustic wave incident from 
0 = 0 . 
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Comparison of Figure 5 to F igure 1 0. 1 0 of Bowman [Ref. 1 6] 
shows a match in both 0 and 5 . This graphical match verifies the 
functioning and output of the differential scattering cross section 
program of Appendix B. 




C = 343 M/s. 

F = 5000 HZ, 

A = .0I0918029O9GI M. 

K«A = I . 

NORMALIZED BCKSCTR = .880073464081 

Figure 5, Sheet 1 of 3; Differential Scattering Cross Section 
Program Output for Comparison to Figure 10. 10 of Bowman [Ref. 16]. 
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CTPG MODULUS 




C = 343 M/S. 

F = 5000 HZ. 

A = .054590]454805 M. 

K*A = 5 . 

NORMALIZED BCKSCTR = 1 . 08 1 2038B888 



Figure 5, Sheet 2 of 3: Differential Scattering Cross Section 
Program Output for Comparison to Figure 1 0. 1 0 of Bowman [Ref. 1 6], 
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5CTRG MODULUS 



10L1 




C = 343 M/S. 

F = 5000 HZ. 

A = . 109180290961 M. 

K*A = 10 . 

NORMALIZED BCKSCTR = .957615718739 



Figure 5, Sheet 3 of 3; Differential Scattering Cross Section 
Program Output for Comparison to Figure 10. 10 of Bowman [Ref. 16], 
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2. Backscattered Cross Section Software Development 



The differential scattering cross section program of Appendix B 
was modified to yield the backscattered cross section program of 
Appendix C. With the modification the normalized differential scattering 
cross section is calculated for the backscattered direction only. Equation 
(86b) is used for this calculation. 

Additionally, the graphics subprograms of the differential 
scattering cross section program were replaced by one graphics 
subprogram. The new graphics subprogram was designed to provide a 
rectangular graph of the normalized backscattered cross section, 

is plotted as a function of ka for a specifiable range of values of ka. 

Figure 6 was generated by the backscattered cross section 
program. It was used for comparison to Figure 10. 1 1 of Bowman [Ref. 16]. 
This comparison served as a verification of the calculated, normalized 
backscattered cross section values. Note the ka axis is labeled K^A. 

Figure 7 was generated to verify the normalized backscattered 
cross section profile was approaching a value of one as ka left the Mie 
scattering regime and entered the "geometric optics" region. Note the ka 
axis is labeled K^A. 
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NORM'D BCKSCTRG CROSS SEIC'N 




Figure 6: Backscattered Cross Section Program Output 
for Verification of Calculated Normalized Backscattered 
Cross Section Values by Comparison to Figure 10. 1 1 of 
Bowman [Ref. 16]. 
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NORM D BCKSCTRG CROSS SEC 




Figure 7; Backscattered Cross Section Program Output 
for Verification of the Normalized Backscattered Cross 
Section Asymptotic Behavior with Increasing ka. 
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3. Acoustic Array Efficiency Software Development 



The acoustic array efficiency program in Appendix D determined 
In general, each major section of Chapter II of this thesis and each 

set of analyzed data is represented by a subprogram of the acoustic array 
efficiency program. 

The acoustic array efficiency program of Appendix D includes the 
non-graphics subprograms present in the backscattered cross section 
program of Appendix C. It also includes subprograms designed to furnish 
stored data, calculate attenuation and range using equations (27) through 
(42), and use developed information to calculate E^^G^ with equation (17). 

Equation (86b) was used to calculate (Tj^. 

The appropriate subprogram and main program calculational 
results were then used in equation (26) to determine E^,. was used in 

equation (92) to determine E^ E^. 

B. DATA ACQUISITION AND SUPPORTING EQUIPMENT 

Data collection occured 28 November 1987, 30 November 1987, and 
7 December 1987. 
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Data collected 28 November and 30 November provided information 



needed for calculation of (r.E E,G„. Data collected 7 December allowed 

D r L U 

calculation of E,G„ and the calculation of E 's. 

t u r 

Collected data is presented in Appendix E. 

1 . Two Way Propagation Path Measurements 

Two way propagation path data were collected 28 and 30 
November. Measurements were used to calculate the product g-j^E^Ej^Gq. 

Figure 8 provides a connection schematic for the data acquisiton 
equipment. 




Figure 8; Data Acquisition Equipment Set-up 
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Four target spheres were available and measured for size. They 
included one hollow aluminum sphere (sphere number one) and three solid 
aluminum and brass spheres. Table ( 1 ) lists the data and the calculated 
values of = Tta^, ka, and 0 -,^. was calculated from equation (86b). 

Table 1: Target Sphere Dimensional Data and Normalized 
Backscattered Cross Section Calculational Results for 
Frequency of 5000 [Hz]. 



Target 

No. 


Sphere 
Diameter [m] 


Sphere 
Radius [m] 


Im^l 


ka 




1 


.2546 ±.0006 


.1273 


.05091 


1 1.62 


.9142 


2 


.1013 ±.0001 


.05065 


.008060 


4.619 


1.000 


3 


,07634±.00003 


.03817 


.004577 


3.479 


1.051 


4 


.06240±.00008 


.03120 


.003058 


2.843 


.6836 



Figure 9 shows the target spheres' theoretical normalized 
backscattered cross section distribution plotted on a graph generated by 
the backscattered cross section program of Appendix C. 
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® 2 4 E 0"" 10 i 2 






Figure 9: Target Spheres' Normalized Backscattered Cross 
Section Distribution as Calculated from Measured Diameters, 
Using an Assumption of Perfect Acoustic Hardness. 



The target spheres were suspended in turn in the anechoic 
chamber. Placement was made using the criteria in Section D. 1 of 
Chapter II. The range R from the array to a target was determined to be 
5.30 [ml 

Temperature T^ and relative humidity Rh along the propagation 
path were determined with a Weathermeasure Temperature and Relative 
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Humidity Meter. was usually about 2 1 [°c] with a relative humidity of 
approximately 50% Atmospheric pressure P was read from a barogram 

supplied by the Naval Postgraduate School meteorological station. 

Specific T^, Rh. and P data are furnished in Appendix E. 

The unshrouded, close packed, hexagonal acoustic array of 
Moxcey [Ref. 8] was used to ensonify each target. 

A = Ng TT rg2 (96) 

determined the aperture area A of the array. Ng = 1 9 was the number of 
speakers in the array and r^ = .0381 [m] was the average speaker horn 
radius [Ref. 8]. 

Moxcey [Ref. 8] determined the half width of the array’s main 
lobe to be 12®. The angular half widths of targets two, three and four 
were all less than 5% of the main lobe's half width. Target one's angular 
half width was 1 1.5% of the main lobe's half width. 

A sinusoidal input from a Hewlett Packard 3314 A function 
generator drove the array via a pre-amplifier. Signal input from the 
HP33 1 4A function generator to the pre-amplifier was set for 20 cycles at 
5000 [Hz]. The voltage input to the pre-amplifier’s power amplifier 
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terminal (see Figure 8) was 3.5 [V] for spheres two, three and four and 3.0 
(V] for sphere one. 

Twenty cycles is considerably less than the number of cycles 
expected in actual echosounder application. It was necessary to keep the 
pulse packet this small because the anechoic chamber’s corner was only 
about 1.5 [m] from the target mounting location. Voltage input for sphere 
one was reduced to 3.0 [V] to avoid clipping of its greater return signal 
voltage. 

Actual transmission voltage inputs to the array from the 

pre-amplifier were measured with a Nicolet 3091 oscilloscope at the 
pre-amplifier’s array terminal. The array was disconnected from the 
pre-amplifier for measurements. values are provided in Appendix E. 

They are also stored in the Trans30 and Trans35 subprograms of the 
acoustic array efficiency program in Appendix D. 

The Nicolet oscilloscope was connected to a Hewlett Packard 
7090 A measurement plotting system. This connection provided a hard 
copy plot of voltage vs. time traces. 

The returned signal was amplified and filtered. 
Pre-amplifier gain Gg was determined asl 1 ,094 by Moxcey [Ref. 8]. 
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Filtering was done with a Wavetek Brickwall filter model 753 A to remove 
the low frequency normal modes of oscillation of the speaker drivers from 
the return signal. Filtering was accomplished at 0 [dB] gain and 5000 [Hz] 
for both low and high pass settings. 

The processed return signal voltages V^, were displayed and 

measured on the Nicolet 3091 oscilloscope. measurements were made 

at the Out 2 terminal of the Wavetek Brickwall filter (see Figure 8). They 
were taken at 10 [p5] intervals for two full cycles along the asymptotic 
region of the return signal trace. 

The Nicolet 3091 oscilloscope was connected to the Hewlett 
Packard 7090 A measurement plotting system. Figure 10 shows the 
Hewlett Packard plotters traces for target sphere two’s return signal 
voltages as printed and subsequently labeled. The figure is representative 
of the oscilloscope traces for each target return. Data for all targets are 
provided in Appendix E. 

The time of return t was determined as the time difference 

r 

between the first positive maxima of the transmitted pulse packet and the 
first positive maxima of the target’s return pulse packet. 
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V i 



V A 



0[V] 




400 [u5] 2,65 [mS] 

2 26 [mS] 

'.8: 27: 'i? 28 NGV S: 



L 

r 



r 




30 95 [mS] 



3452 [mS] 



34 1 3 [mS] 



Figure 10, Sheet 1 of 2; Target Sphere 
Two’s 5000 [Hz] Return Signal Trace 




:8:28:37 28 NOV 87 

Figure 10, Sheet 2 of 2; Target Sphere Two's 
5000 [Hz] Return Signal Trace, Scaled Down 
For Oscilloscope Measurement. 
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Noise voltage measurements were made with the same 
measuring equipment setup used for measurements. The target spheres 
were removed. V„ measurements were taken at a time of return t 

n r 

comparable to the time of measurement of V^. 

measurements were actually made on occasions well 
separated chronologically from the measurements. Therefore, they 

were proportionally corrected in line 2530 of the acoustic array 
efficiency program. The corrections account for the slight variations in 
input voltages and propagation path attenuations that occurred between 
the two measurement occasions. 

Figure 1 1 shows a representative noise trace. Data from the 
noise measurements is presented in Appendix E and stored in the Noise30 
and Noise35 subprograms of the acoustic array efficiency program in 
Appendix D. 

Unless otherwise specified above the remaining data from the 
two way propagation path measurements are stored in the subprograms 
titled Sphere!, 5phere2, Sphere3 and 5phere4. All remaining data are also 
presented in Appendix E. 
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V 




figure 1 1: Representative 5000 [Hz] Noise Trace. 

3.0 [V] Supplied by the HP33 1 4A Function Generator. 

2. One Way Propagation Path Measurements 

One way propagation path data were collected 7 December. 
Measurements were used to calculate with equation (17). 

Data from the one way propagation path measurements are 
provided in Appendix E. Data are also stored in the Gain30 and Gain35 
subprograms of the acoustic array efficiency program of Appendix D. 
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For one way path measurements, i.e. for the determination of 
Ej^Gq, an I vie Electronics Inc. IE 30 A Audio Analyzer (Serial -^805 A 426) 

was used to measure received target intensity 1^ in [dB], Intensity 

measurements were referenced to 10“’^ [W/m^]. The acoustic signal 
supplied for this measurement was a continuous wave at 5000 [Hz]. 
Measurements were taken at the suspension site of the calibrating targets. 

A Fluke 8060 A True RMS Multimeter was used to measure the 
root-mean-square voltage supplied to the acoustic array at its 

input. was squared to determine (V^g)^. 

The Weathermeasure temperature and relative humidity meter 
was unavailable for E^G^ data collection. Hygrothermographs from the 

Naval Postgraduate School's meteorological station were used to estimate 
the anechoic chamber's relative humidity Rh and temperature T^ on 

7 December. Estimates were based on comparisons to the 
hygrothermographs and measured Rh and T^ of 28 and 30 November. 
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C. RESULTS 



Data analysis was performed with the aid of the acoustic array 
efficiency program of Appendix D. Data used in calculations is furnished 
in Appendix E. 

The product generated from two way popagation path 

data and application of equation (26). The average of calculated ctj^'s 
compared to the average of the experimental product Cj^Ej^E^G^ gives a 
normalizing product E^Ej^G^^. The experimental product c-j^E^E^G^ and the 
normalizing product E^E,.G^ are used to calculate experimental o-j^'s. 
These results are supplied in Table 2. 



Table 2; Target Sphere Experimental, Normalized 
Backscattered Cross Section Calculational Results. 



T arget 
No. 


ka [m^] 


Calc'd 


brio 


r t 0 


Exp1 


1 


1 1.62 


.9142 


55.90 




.947 


2 


4619 


1.000 


59.88 




l.Ol 


3 


3.479 


1.051 


61.62 




1. 04 


4 


2.843 


.6836 


37.98 




.644 


Avg. 




.912 


53.85 


59.0 
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NORM'D BCKSCTRG CROSS SEC'N 



The experimental are plotted on a graph generated by the 

backscattered cross section program of Appendix C. This plot is Figure 
12. It is similar to Figure 9 and gives as a function of ka. Note ka is 

labeled K^A. 




° 6^ S -i E e '■ 10 12 



K*fl 



Figure 12; Target Spheres’ Experimental. Normalized 
Backscattered Cross Section Distribution Compared to 
a Calculated Backscattered Cross Section Distribution 
Curve Generated by the Backscattered Cross Section 
Program of Appendix C. 
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As ka increases the experimental distribution rises from just 
below the calculated curve to Just above it, 

was determined using one way propagation path data and equation 
(17). Ej^Gq, or E^, was assumed to be the same for spheres two, three, and 

four because the supply voltage and target locations were the same. The 
angular half widths of target spheres two, three and four were between 
4.6% and 2.8% of the echosounder's main lobe's half width. 

Table 3: Acoustic Array Product E^G^ 

Determined from One Way Propagation 
Path Measurements and Equation (17). 



Target 

No. 


HP3314A 

Func'n Gen'r 
Supply Volt. [V] 




1 


3.0 


94.70 


2 


3.5 


95.1 1 


3 


3.5 


95.1 1 


4 


3.5 


95.1 1 



Calculated E^G^ and two way propagation path data were used in 
equation (26) to determine E^. values were calculated from the target 



60 



spheres' measured diameters, 2a. was calculated as Tia^. Calculation 
results are furnished in Table 4. 

Table 4: E^, and E^ Determined From Two Way Propagation 

Path Data and Equation (26). Ej^E^ Determined From Assumption 
of Efficiency Reciprocity. 



Target 

No. 


V[m^] 


Calc'd 


ffb [m2] 






1 


.05091 


.9142 


.04654 


.6456 


.4169 


2 


.008060 


1.000 


.008060 


.6296 


.3964 


3 


.004577 


1.051 


.004810 


.6167 


.3803 


4 


.003058 


.6836 


.002090 


.5841 


.341 1 



The E^'s determined for Moxcey's unshrouded hexagonal array are in 

the range of .58 to .65. This is above the value range of 0.5 + 0. 1 
determined by Weingartner for his square acoustic array [Ref. 6]. 

Table 4 results were used to construct Figures 1 3 and 1 4. Figure 1 3 
gives E^ as a function of the calculated, normalized backscattered cross 

section, Figure 14 gives E^ as a function of the apparent target size. 



61 



4 



Target Number 







Figure 13; Acoustical to Electrical Power Conversion Efficiencies, 
E as a Function of the Normalized Backscattered Cross Sections, 

(Tjj, Calculated from Measured Sphere Sizes. 



Target Number 



1 




Figure 14; Acoustical to Electrical Power Conversion Efficiencies, 
E^ as a Function of the Apparent Target Sizes, 
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Recalling equation (20), 





G(Q)] [-^] 



(20) 



and equation (12), 



p z 

4jiR 



-CtR . 

• r-. 2 ^ f 



( 12 ) 



serves as a reminder that the received acoustic power, P^g. depends on the 



apparent target size, 

Figure 14s trend and the dependency of P^g on were used as 
suggestions for the construction of Figure 15. Figure 15 is a plot of 



Target Number 

4 3 2 1 




Figure 15; Acoustical to Electrical Power Conversion Efficiencies, 
E^ as a Function of the Received Acoustic Powers, P^g. 
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versus P 



Ppg was calculated using the definition (from equation (21 )) 



ra 



Pra - Pp / Ep 



( 97 ) 



and the relationship (from equation (15)) 



Pr “ / 2 



( 98 ) 



Note Figure 15's trend is for to increase with an increase in P^^ for 



the range of target sizes investigated. 
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IV. CONCLUSIONS AND RFCOnflENDATIONS 



A CONCLUSIONS 

The larger magnitudes of determined for Moxcey's hexagonal array 

suggest it may be more efficient than Weingartner's square array. 

However, the drivers for each array are identical. 

The 0.08 to 0.15 differences in efficiencies are more likely the result 
of the difference in calibration methods. Weingartner's measurements 
[Ref. 6] used steady state, continuous waves. The use of steady state, 
continuous waves allows the development of standing waves. This 
calibration investigation used pulsed waves. The use of pulsed waves did 
not provide sufficient time for standing waves to completely develop in 
the backscattered energy. [Ref. 26] This calibration method should provide 
more accurate efficiency values since it more nearly approximates the 
echosounder system's actual field operating configuration. 

Additionally, steel braid and nylon lines were used to suspend the 
target spheres. These supports probably reflected enough acoustic energy 
to have a small influence on results. Other fixtures in the anechoic 
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chamber, such as support clamps and lights, may have also contributed 
small amounts of reflected energy to the returned signal. 

Calibration of the computer controlled echosounder required 
determination of the product E^E^. The calibration process assumed 

E^ = The calibration results in Tables 3 and 4 for spheres two, three, 
and four indicate E^ varied with constant E^G^, or constant Ej^. Such a 
variation negates the assumption of reciprocity for E^ and Ej^. 

Figure 15 suggests the acoustic array's efficiency E^ increased in 
response to an increase in P and (V ) . 

^ ra rms r 

The acoustic gain for spheres two, three and four was probably not 
constant, as assumed. However, Gq is the maximum for G(0), and E^ is 

inversely proportional to G^, per equation (26). An averaging process 

accounting for the slight differences in gain encountered by spheres two, 
three and four would have increased the data points spread observed in 
Figures 1 4 and 15. 

Additionally, the pre-amplifier's gain may not have been constant for 
the range of received voltage values encountered. This would have 
introduced either an apparent increase in E^ or dampened an increase in E^ 
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with an increase in Such a gain variation could have effectively 

negated the assumption of efficiency reciprocity. However, this 
pre-amplifier is approximately 0.1% linear for less than 10 [V] peak 
[Ref. 26]. Received voltages in this investigation were all less than 10 [Vj. 

B. RECOMMENDATIONS 

Because efficiency reciprocity is now in question it is necessary to 
determine Gq by calculation or measurement. Such a determination would 

allow the desired calibration quantity E^Ej^ to be written as 

(99) 

One method of calculating Gq would involve integrating the acoustic 

array’s calculated or measured intensity pattern (see Moxcey [Ref. 8]) for 
determination of an equivalent isotropic emitter. To determine Gq the 

centroid intensity of the array’s emitted intensity pattern would be 
compared to that of the equivalent isotropic emitter. This would be 
repeated for each transmission power level of interest. 
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Calculation of Gq would allow a determination of using E^Gq . Such 
a determination of E^ would allow a reciprocity comparison to E^ and 
support a more accurate evaluation of E^E^. 

The accuracy of and E^G^ could be further refined as follows. 

Use a calibrating site that allows the use of a pulse train closer in 
length to actual field use. 

Make target return and noise voltages and site intensity 
measurements as chronologically close as possible. 

Measure the barometric pressure in the propagation path directly 
along with the relative humidity and temperature. 

The received and transmitted electrical powers for E^ and the 

transmitted electrical power for E^Gq were all measured at the 

pre-amplifier end of the pre-amplifier/array transmission line. This 
included transmission line power losses due to impedance as efficiency 
decreases. The present calibration process includes the transmission line 
impedance as part of the array impedance. 

The amount this impedance power loss actually affects E^ could be 

evaluated by determining the impedances seen in each direction at the 



68 



pre-amplifier's array terminal. This would be an opportune time to 
verify the validity of transmission line impedance reciprocity. 

Use the Nicolet oscilloscope to record a triggered readout from a 
calibrated microphone suspended at the target site at 45° to the centroid's 
axis. This would allow E^Gq to be determined for a pulse packet rather 

than for a continuous wave. Coupled with a proper calibrating site 
selection this technique should minimize interference problems such as 
anechoic chamber corner reflections. 

A larger number of apparent target sizes spanning the desired range 
of apparent target sizes should be used to redetermine for Moxcey's 

array. The E^, distribution should be examined for a smooth and asymptotic 
behavior. Comparative E^ values for Moxcey's array could be determined 

using the method outlined by Weingartner (Ref. 6]. 

Weingartner's square array's E^ could also be determined with this 

thesis' calibration process and compared to Weingartner's value [Ref. 6). 

The refined E^ curve in the now calibrated echosounder system should 
be used to develop an absolute C.j.2 plot of an air mass. An established 
method such as a tower mounted vertical array of thermocouples would be 
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used simultaneously to generate a comparative absolute plot of the 
same air mass. 

C. SUMMARY 

Realizing the potential of the computer controlled echosounder for 
analyzing lower atmospheric turbulence required a calibration of the 
echosounder. This thesis described theory and software for performing 
calculations crucial to a backscattered cross section calibration of the 
echosounder. The backscattered cross section calibration was evaluated 
in this thesis. 

E^‘s values indicate the evaluated calibration process requires some 

refinement. They also indicate the calibration method possesses 
sufficient merit to warrant further development and that the software 
performs as intended. 
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APPENDIX A 



ANECHOIC CHAMBER DESCRIPTION 

The anechoic chamber is located in room 019 of building 232 at the 
Naval Postgraduate School. Acoustic research on sound sources, sound 
receivers and sound scatterers may be conducted in the chamber \with a 
minimum of interference from wail reflection and external noise. 

Wall reflection is minimized by 102 centimeter deep wedges made of 
P.F. 612 fiberglass. The wedges are attached to the walls, ceiling and 
floor of the chamber. Approximately M2 cubic meters of fiberglass is 
used to absorb 99% of incident sound energy for frequencies greater than 
100 Hertz. External noise isolation is accomplished by separating the 
chamber's inner concrete block sides from the outer 12 inch concrete 
walls with a two inch thick lining of fiberglass and cork. 

A usable region of approximately 8.2 meters by 4.3 meters by 3.4 
meters is available. This region is floored by a grid of 225 wire cables, 
with each cable placed in a tension of 150 to 200 pounds per square inch. ^ 



^This Appendix was compiled from the description posted by the 
entrance to the anechoic chamber instrumentation and control room. 
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APPENDIX B 



DIFFERENTIAL SCATTERING CROSS SECTION PROGRAM 



Examples of the output of this program are Figures 4 and 5, located in 
Chapter III. 



]0 


! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 


20 


! 














30 


REM 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


40 


REM 


ACOUSTICALLY 


HARD 


SPHERE 


SCATTERING OF PLANE ACOUSTIC WAUES 


50 


REM 














60 


REM 


THIS PROGRAM 


CALCULATES 


AND 


GRAPHS THE DIFFERENTIAL SCATTERING 


70 


REM 


CROSS SECTION AND/OR 


SCATTERING 


MODULUS FOR AN ACOUSTICALLY HARD 


80 


REM 


SPHERE OF 


DIAMETER D 


ENSONIFIED 


BY AN INCIDENT PLANE WAUE OF 


90 


REM 


FREQUENCY 


F. 










100 


REM 














130 


REM 


+++++++++++++++++++++++++++++++++ +++++++++++++++++++++++++++++++-F+ 


140 


1 














150 


I++++ UARIABLE 


DECLARATIONS 


■ AND DEFINITIONS ++++++++++++++++++++++++++++ 


160 












j 




170 












j 




180 


REAL 


C 








! 


SPEED OF SOUND. 


1 90 












j 




200 
7 1 (A 


REAL 


Freq 








1 

1 


FREQUENCY OF THE EMITTER.' 


L 1 ^ 
220 


REAL 


A 








j 


RADIUS OF THE SPHERE. 


230 












1 




240 


REAL 


L 








! 


ORDER OF THE FUNCTION CALCULATED 


250 












1 


IN A SUBROUTINE. IT IS USED 


260 












! 


AS A LOOP INDEX AND TO DESIGNATE 


270 












j 


THE ARRAY ELEMENT FOR THE FUNCTION 


280 












1 


OF CORRESPONDING ORDER. 


290 












! 




300 


REAL 


Lnax 








1 


MAX. SIGNIFICANT ORDER OF THE SUM 


310 












! 


DETERMINING THE DIFFERENTIAL 


320 












! 


SCATTERING CROSS SECTION (DSCS) 


330 












! 


AND SCATTERING MODULUS < SM ) . 


340 












! 




350 


REAL 


K 








! 


THE WAUE NUMBER. 2*PI*FREQ/C. 


360 












1 




370 


REAL 


Ka 








j 


K*A. THE ARGUMENT OF THE REGULAR 


380 












! 


AND IRREGULAR SPHERICAL BESSEL 


390 












1 


FUNCTIONS AND THEIR DERIUATIUES. 


400 












! 
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410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 



REAL J(51) I ARRAY FOR THE REGULAR SPHERICAL 

! BESSEL FUNCTIONS ( RSB ) OF 
! ARGUMENT KA AND ORDER 0 THROUGH 
I LMAX. 



REAL Dj(5I ) 



ARRAY FOR THE DERIUATIUES OF THE 
RSB (THE DRSB). 



REAL Y(51 ) 



ARRAY FOR THE IRREGULAR SPHERICAL 
BESSEL FUNCTIONS (ISB) OF 
ARGUMENT KA AND ORDER 0 THROUGH 
LMAX. 



REAL Dy(51 ) 



ARRAY FOR THE DERIUATIUES OF THE 
ISB (THE DISB). 



REAL P(51 ) 



ARRAY FOR THE LEGENDRE 
POLYNOMIALS OF ARGUMENT COS(PHI). 



REAL Sn(360) 



THE ARRAY REPRESENTING THE 
SCATTERING MODULUS ( SM ) . 



REAL D5(360) 



THE ARRAY REPRESENTING THE 
DIFFERENTIAL -SCATTERING CROSS 
SECTION NORMALIZED BY THE CROSS 
SECTIONAL AREA OF THE SPHERE 
(DSCS). 

NOTE DS(PHI ) = 4*SM(PHI )"2/KA''2. 



++++ NOTE ++++ 

THE ARRAY DIMENSIONS (30+21) NEED 
TO BE INCREASED TO ACCOMODATE 
KA>30. 



I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



I++++ main program +++++++++++++++++++++++++++++++++++++++++++++++++++++ 



CALL 

CALL 

CALL 



Init(C.Freq.A.Lnax ,Ka) ! INPUT AND CALCULATE REQUIRED 

I PARAMETERS. 

! 

Di fsctr( Lnax ,Ka ,Dj( ♦ ) ,J( ♦ ) ,Dy( • ) ,Y( * ) ,P( * ) ,Sn( * ) ,Ds( ♦ ) ) 

! CALCULATE THE SM( ♦ ) AND 0S(O. 



Outpt(C,Freq ,A,Ka,Sn( ♦ ) ,Ds( • ) ) 



OUTPUT THE SM(») AND/OR DS(M. 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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910 

320 

330 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 



MAIN PROGRAM CONCLUSION ++++++++++++++++++++++++++++++++++++++++++ 

! 

PRINTER IS 1 
PRINT 

PRINT "END*' I ADUISES THE USER OF THE MAIN 

PRINT ! PROGRAM’S CONCLUSION. 

END 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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1030 

1040 

1050 

t 060 

1070 

1080 

1090 

1 100 

1110 

1 120 

1 130 

1 1 40 

1 150 

1 160 

1170 

1 180 

1 1 90 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1 440 

1450 

1 460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Ini t(C ,Freq ,A .Lmax ,Ka ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE REQUESTS THE APPROPRIATE INPUT VARIABLES TO DETERMINE 
REM C, FREQ, A, LMAX , AND KA. 

REM ++++ + + +++ + + ++++++ + +++++ + +++++ + + +++ + ++++++ + + +++ + +++-»-++ + + + +++ + + + ++ + + 

! 

!++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL Temp 
REAL Circ 
REAL Diam 



AMBIENT AIR TEMPERATURE IN DEG C. 
CIRCUMFERENCE OF THE TGT. SPHERE. 
DIAMETER OF THE TARGET SPHERE. 



! + + +++ +++++++ + + +++-I-+ + + + + + + + H 



h+++++++++++++++++++++++++++H 



!++++ REQUESTING AND CALCULATING SYSTEM PARAMETERS ++++++++++++++++++++++ 



PRINTER IS 1 
Choice=0 

WHILE ChoiceOl AND Choice<>2 
PRINT 

PRINT ‘'DO YOU DESIRE.TO INPUT" 
PRINT " (1 ) THE SPEED OF " ; 

PRINT "SOUND," 

PRINT " OR (2) THE AMBIENT AIR " 
PRINT "TEMPERATURE?" 

PRINT 

INPUT Choice 
END WHILE 

IF Choice=1 THEN 

PRINT "PLEASE ENTER C IN M/S." 
INPUT C 
END IF 

IF Choice=2 THEN 

PRINT "PLEASE ENTER TEMP IN "? 
PRINT "DEGREES CELSIUS. " 

INPUT Temp 
Temp=Temp+273 . 15 
C=20.05«SQR(Temp ) 

END IF 

PRINT 

PRINT "PLEASE ENTER THE ECHOSOUNDER 
PRINT " FREQUENCY IN HZ." 

PRINT 

INPUT Freq 



INPUTS ARE REQUESTED ON THE CRT. 
DETERMINING C 



! DIRECTLY ENTER C, 



! OR CALCULATE C FROM TEMP. 



! *** DETERMINING FREQ 
! DIRECTLY ENTER FREQ. 
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« « « 



1580 

1590 

1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1 800 
1810 
1820 
1830 
1 840 
1850 
1860 
1870 
1880 
1890 
1 900 
1910 
1920 
1 930 
1 940 
1950 
! 960 
1 970 
1 980 
1 990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 



Choice=0 ! DETERMINING KA 

WHILE ChoiceOl AND Choice<>2 AND Choice<>3 
PRINT 

PRINT "DO YOU DESIRE TO INPUT " 

PRINT " ( 1 )THE SPHERICAL " ; 

PRINT "BESSEL ARGUMENT K*A." 

PRINT " (2 ) THE SPHERE’S " i 

PRINT "CIRCUMFERENCE," 

PRINT " OR (3) THE SPHERE’S 
PRINT "DIAMETER?" 

PRINT 

INPUT Choice 
END WHILE 

1 

IF Choice-1 THEN ! DIRECTLY ENTER KA , 

PRINT "PLEASE ENTER K*A. " 

INPUT Ka 
END IF 



IF Choice-2 THEN 

PRINT "PLEASE ENTER THE TARGET " 
PRINT "CIRCUMFERENCE IN CM." 
INPUT Circ 

Ka=Circ*( 10^(-2))*Freq/C 
END IF 

IF Choice=3 THEN 

PRINT "PLEASE ENTER THE TARGET " 
PRINT "DIAMETER IN CM." 

INPUT Dian 

Ka-Dian*( 10"(*-2 ) )*PI*Freq/C 
END IF 



K=2*PI*Freq/C 

A=Ka/K 

Lna/- INT ( Ka + 2 1 ) 



! CALCULATE KA FROM THE 
! TARGET CIRCUMFERENCE, 



! OR CALCULATE KA FROM THE 
! TARGET DIAMETER. 



! *** DETERMINING REMAINING *** 

! SYSTEM PARAMETERS 
! CALCULATE THE WAUE NUMBER. 

I CALCULATE THE SPHERE RADIUS. 

! MAX. SUMMATION ORDER CALCULATED 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

I++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 



! 

+ + + + i + + + + + + + + + + + + + + + 4 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2550 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 



) $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Difsctr{LnaK,Ka,Dj<* ),J( *),Dy(* >,Y(* ),P( * ) ,Sn(* >,Ds(* >) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE ARRAYS SM( * > AND DS( * > REPRESENTING THE 

REM SCATTERING MODULUS AND THE DIFFERENTIAL SCATTERING CROSS SECTION 

REM NORMALIZED BY THE CROSS SECTIONAL AREA OF THE SPHERE, RESPECT lUELY . 
REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

} 

!++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



FACTOR COMMON TO RSM AND ISM. 

THE REAL COMPONENT OF THE 
SCATTERING MODULUS, SM. 

THE IMAGINARY COMPONENT OF THE 
SM. 

FACTOR COMMON TO SM(PHI) & DS(PHI). 

THE (POLAR) ANGLE OFF THE AXIS 
OF PROPAGATION OF THE INCIDENT 
PLANE UAUE, WITH ORIGIN AT THE 
SPHERE CENTER. COS(PHI) IS THE 
ARGUMENT OF THE LEGENDRE 
POLYNOMIALS P( ♦ ) OF ORDER 0 
THROUGH LMAX. 






REAL Fac 
REAL Rsn 

REAL Ism 

REAL Fact 
REAL Phi 



!++++ CALCULATING THE SM AND DSCS ARRAYS +++++++++++++++++++++++++++++++ 



CALL Drsb ( Lmax ,Ka ,D j < ♦ ) , J ( ♦ > > 

CALL D 1 sb ( Lmax ,Ka ,Dy < * > , Y( ♦ > > 

FOR Phl=0 TO 360 

CALL Leg(Lmax ,Ka ,Phi ,P( * ) ) 

R5m=0 
I 5m=*0 

FOR L=0 TO Lnax 



! CALCULATE THE DERIUATIUES OF THE 
! REGULAR SPHERICAL BESSEL FUNCTIONS 
! OF ORDER 0 THROUGH LMAX. 

! CALCULATE THE DERIUATIUES OF THE 
! IRREGULAR SPHERICAL BESSEL FUNC^NS 
! OF ORDER 0 THROUGH LMAX. 



CALCULATE THE ARRAY P(*). 

CALCULATING THE REAL AND IMAGINARY 
COMPONENTS OF THE SM. 



Fac = (Dj(L )M2*L+1 >*P( L ) >/( Dj ( L >*Dj ( L > + Dy( L >*Dy( L > ) 
R5P'i = Rsm+(Fac*Dy(L > > 
l5m-Ism + (FdC*Dj(L > > 

NEXT L 
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2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2850 



Fact=R5n*R5nrMspn*l5n 



Sn(Phi ) = 50R(Fact ) 

Ds(Phi ) = 4*(Fact )/(Ka*Ka ) 

NEXT Phi 



CALCULATING THE ARRAY REPRESENTING 
THE SM(PHI ). 

CALCULATING THE ARRAY REPRESENTING 
THE DIFFERENTIAL SCATTERING CROSS 
SECTION NORMALIZED BY THE CROSS 
SECTIONAL AREA OF THE SPHERE. 



I 

I 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

j 

SUBEND 



I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
! 
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2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

31 10 

3120 

3130 

3140 

3150 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Drsb ( Lmax , Ka , D j ( * ) , J ( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE REGULAR SPHERICAL 
REM BESSEL FUNCTIONS (DRSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . 

REM THE CALCULATION OF THE DERIUATIUE INUOLUES THE RSB’S OF THE SAME 

REM AND PREVIOUS ORDER. 

REM ++++++++++++++ +++++++++++++++++++++++++++++++++++++++++++++ +++++++ 
! 

!++++ CALCULATING THE DERIUATIUES OF THE RSB +++++++++++++++++++++++++++ 

I 



CALL Rsb( Lnax ,Ka , J ( ♦ ) ) 

RAD 

Dj(0 )-(“J(0 )+COS(Ka ) )/Ka 
FOR L=1 TO Lnax 

Dj(L ) = "( (L+1 )/Ka )*J(L) + J(L-1 ) 
NEXT L 



I CALCULATE THE RSB UALUES J(0. 

! ANGLE UNITS FOR DRSB CALC’NS, 

! CALCULATE THE INITIAL DRSB DJ(0). 
! CALCULATE THE REMAINING DRSB. 

I 

I 



i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



! 

!+ + + + MODULE CONCLUSION + + + + + + + + + + + + + + + + + + + + + + +++ + + +++++ + +++ ++++ + -I- + + + + + + 

I 

SUBEND 



1 



!+ + + +++ + + + + + + + + + + + + +++ + +++ + + + + + + +++ +++++ + + + +++ + + + + +++ + + + + + + + + + + + + + + + + + -I- + 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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31G0 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

32G0 

3270 

3280 

3290 

3300 

3310 

3320 

3350 

3340 

3350 

3360 

3370 

5380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

34G0 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

35G0 

3570 

3580 

3590 

3600 

3G10 

3G20 

3B30 

3G40 

3G50 

36G0 

3G70 

3G80 

3G90 

3700 



! 

j $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB RsbCLnax ,Ka ,J( ♦ ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE REGULAR SPHERICAL BESSEL FUNCTIONS 
REM (RSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION USES THE CONTINUED FRACTION APPROACH, 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

j 

!++-f+ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL 


R1 


! 

1 

1 THE EUOLUING RATIO OF J(0)/J(L). 
1 


REAL 


Nunr 


! THE EVOLVING NUMERATOR FACTOR USED 

! IN CALCULATION OF RL. 

1 


REAL 


Oenr 


! 

! THE EUOLUING DENOMINATOR FACTOR 
! USED IN CALCULATION OF RL. 

1 


REAL 


Scexp 


! EUOLUING SCALING EXPONENT FOR UERY 
1 LARGE OR UERY SMALL VALUES OF RL. 

1 


REAL 


Nf lag 


! 0 FLAG FOR NUMR. 

j 


REAL 


Of lag 


! 0 FLAG FOR OENR. 


REAL 


Bnm 


! EVOLVING TERM USED IN THE 
! CALCULATION OF NUMR AND OENR. 
1 


REAL 


B i nc 


! INCREMENT USED IN THE EVOLUTION 
! OF BNM. 

1 

! 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ CALCULATING THE RSB ++++++++++++++++++++++++++++++++++++++++++++++ 

j 

j 



RAD 



! ANGLE UNITS FOR CALC’NS. 



J(0) = (SIN(Ka) )/Ka 



FOR L«1 TO Lnax 
Rl = l .0 
Nunr»0 . 

Oenr«0 . 
Scexp"0 , 

Nf lag=0 . 

D f 1 ag=0 . 

Bnn» 1 . 0/Ka 
B 1 nc«2 . 0/Ka 



I 

I *** CALCULATE THE INITIAL *** 
I RSB, J(0). 

! 

! CALCULATING THE RATIO 

I J(0)/J(L). 

! INITIALIZING RL. 

I INITIALIZING NUMR. 

! INITIALIZING OENR. 

I INITIALIZING SCEXP. 

1 INITIALIZING NFLAG. 

! INITIALIZING OFLAG, 

! INITIALIZING BNM. 

! CALCULATE BINC. 
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3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 



FOR Ordcnt=l TO L 
Bnn=Bnn+B 1 nc 
IF Nflag=0 THEN 
Nunr=Bnn“Nunr 
IF Nunr<>0 THEN 
R1=R1 *Nunr 
Nunr= 1 . 0/Nunr 

IF ABS(R1 )>10"250 THEN CALL 
ELSE 

Nf lag=l 
END IF 
ELSE 

Nf lag=0 
END IF 
NEXT Ordcnt 



WHILE NimrOOenr OR Nflag=l OR Of 
Bnn=Bnn+Binc 

IF Nflag=0 THEN 
Nunr=Bnn~Nunr 
IF Nunr<>0 THEN 
R1=R1 *Nunr 
Nunr=^ 1 . 0/Nunr 

IF ABS(R1 )>10"250 THEN CALL 
ELSE 

Nf lag- 1 
END IF 
ELSE 

Nf 1 ag=0 
END IF 

IF Oflag=0 THEN 
Denr=Bnn~Denr 
IF Denr<>0 THEN 
R 1 =R 1 /Denr 
Denr= 1 , 0/Denr 

IF ABS(R1 )>10"25O THEN CALL 
ELSE 

D f 1 ag= 1 
END IF 
ELSE 

Df lag=0 
END IF 

END WHILE 



! * CALCULATE THE L ‘’NAKED** 

! * NUMR FACTORS OF J(0)/J(L). * 

! INCREMENT BNM . 

! CHECK NUMR*0 FLAG NOT SET; 

! UPDATE NUMR IF NUMRO0. 

! UPDATE RL IF UPDATED NUMRO0. 

! PREPARE NUMR FOR NEXT EVOLUTION. 
! CHECK SCALING OF UPDATED RL « 

Sea ie( R 1 ,Scexp ) 



! IF NFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



! ♦ CALCULATE THE REMAINING NUMR * 
! ♦ AND DENR TERMS OF J(0)/J(L). * 
! NOTE WHEN NUMR=OENR THEN 
! RL=J(0)/J(L). 

ag= 1 

! INCREMENT BNM. 

! 

! CHECK NUMR-0 FLAG NOT SET; 

! UPDATE NUMR IF NUMRO0. 

I UPDATE RL IF UPDATED NUMRO0. 

! PREPARE NUMR FOR NEXT EUOLUTION. 
! CHECK SCALING OF UPDATED RL. 

Sea le( R1 ,Seexp > 



I IF NFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



! CHECK DENR-0 FLAG NOT SET; 

I UPDATE DENR IF DENRO0. 

! UPDATE RL IF UPDATED DENRO0. 

! PREPARE DENR FOR NEXT EUOLUTION. 
! CHECK SCALING OF UPDATED RL . 

Sea 1 e ( R 1 .Seexp > 



! IF DFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 
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4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 



J(L)»J(0 )/(Rl*10''(5cexp ) ) ! J(L) DETERMINED. *** 

1 

NEXT L 

I 

! 

! ++++++++++++++++++ 4 -+ +++++++++++++++++++++++++++++++++++++++++++++++++++ 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

I 

I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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4410 

4420 

4430 

4440 

4450 

44B0 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

45G0 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Sea 1 e( R 1 ,Scexp ) 

REM +++++++++++++++++++++++++++++++ 4 -++++++++++++++++++++++++++++++++++ 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN THE 
REM ONGOING CALCULATION OF RL . 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++ SCALING RL +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



! 



IF Rl>10"250 THEN 
Rl=Rl»10''(-250 ) 
Scexp=Scexp+250 
END IF 

IF RK10M-250) THEN 
R1=R1*10"250 
Scexp^Scexp-250 
END IF 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

i 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++H 

! 

SUBEND 



+++++++++++++++++++++++++++ 4 +++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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4720 

4730 

4740 

4750 

47B0 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

48G0 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

49B0 

4970 

4980 

4990 

5000 

5010 

5020 

5030 



! 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 



SUB Di sb ( Lnax ,Ka ,Dy( * ) , Y( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE IRREGULAR SPHERICAL 
REM BESSEL FUNCTIONS OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION OF THE DERIUATIUE USES THE ISB’S OF THE SAME AND 
REM PREVIOUS ORDER. 

REM + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +++ + +++ + + + + + + + + + + -♦■ + + 

! 

!++++ CALCULATING THE DERIUATIUES OF THE I5B +++++++++++++++++++++++++++ 

! 



CALL l5b(Lmax ,Ka ,Y( O ) 

RAD 

Dy(0) = (-Y(0) + SIN(Ka>)/Ka 
FOR L=1 TO Lnax 

Dy(L)=-((L+l )/Ka >*Y(L) + Y(L-1 ) 
NEXT L 



! 

I CALCULATE THE ISB UALUES . 

I 

! ANGLE UNITS FOR CALC’NS. 

! THE INITIAL ISB DERIUATIUE. 

'! THE REMAINING ISB DERIUATIUES. 



! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++ MODULE CONCLUSION ++-M-++++++++++++++++++++++++++++++++++++++++++++ 
SUBEND 

!+ + + + + + + + + + + + + -t- + + + + + + + + + + + + + + + + + + + + + + i- + + + + + i- + + + + + + + + + -M- + + + + + + + + + + + + + + + + + 
! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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5040 

5050 

5060 

5070 

5080 

5090 

5100 

51 10 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Isb(Lndx ,Ka ,Y( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE IRREGULAR SPHERICAL BESSEL FUNCTIONS 
REM OF ARGUMENT KA AND ORDER L THROUGH LMAX . THE CALCULATION IS 
REM ACCOMPLISHED BY FORUARD RECURSION. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



!++++ CALCULATING THE ISB ++++++++++++++++++++++++++++++++++++++++++++++ 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 



RAD 

Y(0 ) = -COS(Ka )/Ka 

Y( 1 )*=(Y(0)-SIN(Ka))/Ka 

FOR L-2 TO Lnax 



ANGLE UNITS FOR CALC’NS. 
INITIAL ISB’S CALCULATED. 



! REMAINING ISB’S CALC*D. 



Y(L)-((2*L-1 )/Ka)*Y<L-l )-Y<L-*2) 
NEXT L 



SUBEND 
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5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

56 10 

5620 

5630 

5640 

5650 



1 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Leg( Lnax ,Ka , Ph 1 ,P ( * ) ) 

REM +++ +++ +++ + + + +++ + + + + + ++++++ + +++-}-+++++++ + + +++ + + + ++++++++++ + + +++ ++++ + 

REM THIS MODULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT 
REM COS(PHI) ANO ORDERS 0 THROUGH LMAX. CALCULATION IS ACCOMPLISHED 
REM BY FORUARO RECURSION. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++ CALCULATING THE LEGENDRE POLYNOMIALS +++++++++++++++++++++++++++++ 

I 



DEG 

X=COS(Phi ) 



ANGLE UNITS FOR CALC’NS. 
ARGUMENT OF THE LEG. POL’S. 



P(0)-1 ! INITIAL LEG. POL’S. CALC’D. 

P( 1 ) = X 

FOR L=2 TO Lnax ! REMAINING LEG. POL’S. CALC’O. 

P(L )-< 2-'( 1 /L ) )*X*P(L-1 )-( 1>< 1 /L ) )»P(L-2 ) 

NEXT L 



! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

1 $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
1 
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5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6 1 30 

6140 

6150 

6160 

6170 

6180 

6190 



! 

j $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB 0utpt(C ,Freq .A.Ka ,Sn< * ) ,Ds( * ) ) 

REM + + +++ + + + + + + + + + +++ + + + + + +++ + + + + +)- + + + + +++ + + ++++ + + + + + + + ++++ +++ + ++ + + + + + 

REM THIS MODULE CONTROLS THE OUTPUT OF THE CALCULATED DSCS AND SM IN 

REM GRAPHICAL PRESENTATIONS. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL Desire 



INTERACTIUE LOOP CONTROL. 



REAL Rpt 



POLARPLT SCALING CONTROL! A 
REPETITION INDEX. 



+++++++++++++++++++++ +++++H 



I- + + + + + + + + + + + +H 



++++ DETERMINING THE DESIRED OUTPUT +++++++++++++ 



++++++++++ 



Des ire=0 

WHILE DesireOl AND Desire<>2 AND Desire<>3 
PRINT "DO YOU DESIRE- 
PRINT " (1) A PLOT OF THE DSCS," 

PRINT " (2 ) A PLOT OF THE SM , " 

PRINT " OR (3) BOTH?" 

INPUT Desire 
END WHILE 
PRINTER IS 701 



WHILE Desire<4 

IF Desire=l OR Desire=3 THEN 
FOR Rpt-1 TO 2 

CALL Polarplt(Ds( * ) ,Rpt ) 
PRINT 
NEXT Rpt 
END IF 

IF Desire=2 THEN 



CALL 


Recp 1 1 ( Sn( * ) ) 


PRINT 






END IF 






PRINT " 


","C =" 


iC;"M/S." 


PRINT " 


" ,"F =" 


;Freq; "HZ. " 


PRINT " 


" . " A = " 


?A? "M. " 


PRINT ■’ 


, " K • A 


' i K a ! " . " 


PRINT " 


“."NORMALIZED 


PRINT " 


BCKSCTR 


= " ?Ds( 180 ) 



PRINT CHR$( 12 ) 

IF Desire=l OR Desire=2 THEN 
IF Desire=3 THEN Desire=2 
END WHILE 



! PLOTTING THE DSCS. 



! PLOTTING THE SM. 



! 

! ECHO SIGNIFICANT CALCULATED UALUES 
! AND INPUT PARAMETERS. 



! PRINTER FORM FEED. 
! 

Des i re=4 



! 
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G200 

G210 

G220 

G230 

G240 

G250 

G2G0 

6270 

6280 

G290 

G300 



1 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
1 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++■+++++++++ 

I 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Polarpi t (Ds( * ) ,Rpt ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE MAPS THE CALCULATED. NORMALIZED DSCS ARRAY DS( O ON 

REM LINEAR POLAR PLOTS. THE FIRST PLOT INCLUDES THE ENTIRE DSCS. 

REM UHILE THE SECOND PLOT SHOWS THE SCATTERING DETAIL AROUND THE 
REM NORMALIZED RANGE. PLOT SELECTION IS GOUERNED BY A REPETITION 
REM INDEX. Rpt. AND IS CONTROLLED IN THE CALLING MODULE (OUTPT). 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL 


Xgunax 


! 

1 

! THE MAXIMUM ABSCISSA UALUE IN 
! GRAPHICS DISPLAY UNITS. 

1 


REAL 


Ygunax 


1 

! THE MAXIMUM ORDINATE UALUE IN 
! GRAPHICS DISPLAY UNITS. 

1 


REAL 


G u rn a X 


! THE MINIMUM OF XGUMAX AND YGUMAX ? 
1 USED TO DETERMINE PLOTTING AREA. 

1 


REAL 


D t nx 


1 THE MAX. UALUE OF THE ARRAY DS( O 
1 


REAL 


Endp t 


! THE MAX. RADIAL SCALE VALUE IN 
! USER DEFINED UNITS. 

) 


REAL 


Ring 


! GRID RING SPACING IN USER UNITS. 
1 


REAL 


Tick 


! RADIAL AXIS TICK MARK SPACINGS IN 
! USER DEFINED UNITS. 

! 



t 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



!++++ INITIALIZING THE PLOTTER +++++++++++++++++++++++++++++++++++++++++ 

! 

GINIT 

PLOTTER IS CRT .“INTERNAL” 

GRAPHICS ON 
GCLEAR 

I 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6980 

6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

71 10 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 



LABELLING +++++++++++++++++++++++++++++++++ 

I 



1 

I++++ GRAPH CONSTRUCTION AND 



Xgunax=100*NAX( 1 .RATIO) 
Ygunax-100*MAX( 1 ,1 /RAT 10) 
6unax=MIN( Xgunax .Ygunax ) 

CSI2E 3 
LINE TYPE 1. 

MOUE Gunax ,Gufnax/2 
LORG 8 
LABEL “0“ 

MOUE Gunax/2 .Gunax 
LORG 6 

LABEL “PI/2" 

MOUE 0 , Gunax/2 
LORG 2 
LABEL “PI" 

MOUE Gunax/2 ,0 
LORG 4 

LABEL “3*PI/2“ 

UIEUPORT 5 ,Gunax“5 ,5 ,Gunax-5 



Dt nax*=0 
IF Rpt=l THEN 
FOR 1=0 TO 360 

IF Dtnx<Ds(I) THEN Dtnx=Ds< I ) 
NEXT I 
END IF 

Endpt = INT( Dtnx+ 1 ) 

IF Endpt<2 THEN Endpt=2 
SHOD -Endpt .Endpt ,-Endpt .Endpt 

LINE TYPE 4 
R i ng= INT( Endpt /5+ 1 ) 

T 1 c k= . 2 *R 1 ng 

AXES Tick .Tick .0 ,0 ,5 .5 ,2 

LINE TYPE 3 
DEG 

LORG 7 

FOR R=1 TO Endpt/Ring 
LINE TYPE 1 
NOUE 0,R*Ring 
LABEL R«Ring 
LINE TYPE 3 
FOR Angle=0 TO 360 
PLOT 

NEXT Angle 
PENUP 
NEXT R 



! LABELLING MAJOR RADIALS ** 



! ** DEFINES THE TOTAL PLOTTING ** 
I ** AREA IN GRAPH UNITS, *♦ 

I SCALING THE GRAPH TO ** 

! USER UNITS. ** 

I DETERMINING MAX. RADIUS UALUE 
! REQ’D FOR THE PLOT. 



I ISOTROPICALLY SCALES THE GRAPH. 
1 

} ** MAJOR RADIALS CONSTRUCTION 



I 

I RANGE RING CONSTRUCTION ** 



I LABELLING THE RANGE. 
I DRAUING THE RING, 



R*Ring«COS( Ang le ) .R«Ring«SIN (Angle ) 



90 



7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7B4C 

7650 

7660 



! 

FOR Angle=30 TO 150 STEP 30 I ** 30 DEGREE RADIALS ** 

IF Angle<>90 THEN ! CONSTRUCTION, 

PLOT Endpt *COS( Angle ) ,Endpt*SIN( Angle ) ,-2 
DRAW Endpt* COS ( 180+Angle ) ,Endpt*SIN( 180 + Angle ) 

END IF 
NEXT Angle 
PENUP 

I 

! 

1 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++ 
1 

!++++ PLOTTING THE CALCULATED DSCS +++++++++++++++++++++++++++++++++++++ 



LINE TYPE 1 
DEG 

FOR 1=0 TO 360 ! PLOTTING RECTANGULAR COORDINATES 

PLOT D5( I )*COS( I ) .Ds( I )*SIN( I ) ! IN POLAR FORMAT. 

NEXT I 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



+ + + + PRINTING THE PLOT + + + + +++ + + + + + + + + + + + + + + + + + + + + +++ + + +++++ + + + + + + + + + -I--I- 



DUMP GRAPHICS #701 



! 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

I++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

GRAPHICS OFF 
SUBEND 

I 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
! 
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7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 

7850 

7840 

7850 

7860 

7870 

7880 

7890 

7900 

7910 

7920 

7950 

7940 

7950 

7960 

7970 

7980 

7990 

8000 

8010 

8020 

8030 

8040 

8050 

8060 

8070 

8080 

8090 

8100 

81 10 

8120 

8130 

8140 

8150 

8160 



SUB Recp 1 1 ( Srn( ♦ ) ) 

REf1 + + + -I- + + + + + + + + + + + + ++ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +++ 

REM THIS MODULE MAPS THE CALCULATED SCATTERING MODULUS ARRAY SMC ♦ ) ON 

REM A SEMI-LOG PLOT. 



REM ++++++++++++++++++++++++++++ 
! 

!++++ MODULE VARIABLE DECLARATIONS 

REAL Xgurnax 

REAL Ygunax 

REAL Dtnx 
REAL Dtnn 
REAL Mxy 

REAL Mny 



+++++++++++++++++++++++++++++++++++++ 
AND DEFINITIONS +++++++++++++++++++++ 



! THE MAXIMUM ABSCISSA UALUE IN 
I GRAPHICS DISPLAY UNITS. 

I 

! THE MAXIMUM ORDINATE UALUE IN 
! GRAPHICS DISPLAY UNITS. 

! THE MAX. UALUE OF THE ARRAY SM(»). 

! THE MIN. UALUE OF THE ARRAY SMCO. 

! THE MAX. ORDINATE SCALE UALUE IN 
! USER DEFINED UNITS, 
i 

! THE MIN. ORDINATE SCALE UALUE IN 
! USER DEFINED UNITS. 



! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++ INITIALIZING THE PLOTTER +++++++++++++++++++++++++++++++++++++++++ 

I 

GINIT 

PLOTTER IS CRT /'INTERNAL*' 

GRAPHICS ON 
GCLEAR 

I 

j +++ + + + -I- + + + + + + + +++++++++ + + + + + ++++ + ++ + + + + + + + + + + + + ++++ + + +++++++ +++++ + + +++ + 

I 

!+ + + + GRAPH CONSTRUCTION AND LABELLING + + + + + + + + +++ + + + + ++++ ++++++ + +++-1- + + + 



Xgunax=l00*MAX( 1 , RAT 10 ) 
Ygumax=100*MAX( 1 J /RATIO ) 

CSIZE 4 

L ORG 4 

MOUE Xgunax/2,0 
LABEL “PHI (DEGREES)’* 



! 



! 

! 

! ** LABELLING THE ABSCISSA AXIS ** 
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8170 

8180 

8190 

8200 

8210 

8220 

8230 

8240 

8250 

82G0 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8380 

8370 

8380 

8390 

8400 

8410 

8420 

8450 

8440 

8450 

8480 

8470 

8480 

8490 

8500 

8510 

8520 

8530 

8540 

8550 

8580 

8570 

8580 

8590 

8800 

8810 

8820 

8830 

8840 

8850 

8880 

8870 

8880 

8890 

8700 

8710 



DEG 

LDIR 90 
LORG 6 

MOlJE 0 ,Ygumax/2 
LABEL "SCTRG MODULUS’' 

LDIR 0 

UIEUPORT 15,Xgunax-5,l0,Ygunax-5 



Dt nx=0 

Dt mn = Sn( 1 80 ) 

FOR 1-0 TO 180 

IF Dtnx<Sn(I) THEN Dtnx=Sn( I ) 
IF Dtnn>Sn(I) THEN Dtnn=Sn( I ) 
NEXT I 

Mxy = INT(LGT(Dtnx ) >+l 
Mny-INT( LGT( Dt nn ) ) 

UINDOU 0 , 180 ,Mny ,Mxy 

AXES 30,1 .0,Mny ,1 ,1,1.5 

LINE TYPE 3 
FOR I=Mny TO Mxy 
FOR J-1 TO 9 

MOUE 0 ,I+LGT( J ) 

PLOT 180 ,I + LGT( J ) ,-l 
PENUP 
NEXT J 
NEXT I 

FOR 1=80 TO 180 STEP 60 
MOUE I ,Mny 
PLOT I ,Mxy ,-l 
PENUP 
NEXT I 

CLIP OFF 
LINE TYPE 1 

LORG 8 

FOR 1=0 TO 180 STEP 30 
MOUE I ,Mny 
LABEL I 
NEXT I 

LORG 8 

FOR I=Mny TO M>. y 
MOUE 0.1 
LABEL lO'^I 
MOUE 0,I+LGT(5) 

LABEL ( 10"( I + l ) )/2 
NEXT I 



! LABELLING THE ORDINATE AXIS ** 



! ** GRAPH CONSTRUCTION 
I DEFINES THE TOTAL PLOTTING AREA 
! IN GRAPH UNITS. 

! ORDINATE USER UNITS ** 

! ** DETERMINED. 



ANISOTROPICALLY SCALES THE GRAPH 
TO USER UNITS. 

PLOTS THE AXES. 



I ♦ PLOTS THE HORIZONTAL 
! * LOGARITHMIC GRID LINES I.A.W. ♦ 
! ♦ THE VERTICAL SCALE. 



) 

! ♦ PLOTS THE VERTICAL GRID LINES • 
1 * I.A.U. THE HORIZONTAL SCALE. * 



) 

1 ** LABELLING THE ABSCISSA SCALE ** 



) ** LABELLING THE ORDINATE SCALE •• 



! POUER OF 10 LABELS. 

! HALF THE NEXT POWER OF 10 LABELS. 
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8720 

8730 

8740 

8750 

8780 

8770 

8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 

8930 

8940 

8950 

8960 

8970 

8980 

8990 

9000 

9010 



! 

CLIP ON 

I 

! 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ PLOTTING THE CALCULATED SCATTERING MODULUS +++++++++++++++++++++++ 



LINE TYPE 1 
FOR 1=0 TO 180 

PLOT I ,LGT(Sm( 1 + 180 > ) 
NEXT I 



! THIS PLOT IS FOR THE PLANE UAUE 
! INCIDENT FROM THE RIGHT (0 OEG’S). 

! USE SM< I ) IN PLACE OF SMC 1 + 180) 

( FOR A PLANE UAUE INCIDENT FROM THE 
! LEFT ( 180 OEG’S >. 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ PRINTING THE PLOT ++++++++++++++++++++++++++++++++++++++++++++++++ 



DUMP GRAPHICS #701 



j 

!++-M'+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



! 

! ++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

GRAPHICS OFF 
SUBEND 



j 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



I 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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APPENDIX C 



BACK5CATTERED CROSS SECTION PROGRAM 



Examples of the output of this program are Figures 6 and 7, located in 
Chapter III. 



!0 


1 ^ ^ ^ ^ ^ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ 


20 


! 






30 


REM 


++ + + + + + + + + + + -I- + + +++ + + + -I- + -I- ++-I- + + + -I- + -I- + + + ++-I- + + + + + + + + +++ +++ + + + + + +++ + + + + + 


40 


REM 


ACOUSTICALLY 


HARD SPHERE BACKSCATTER ING OF PLANE ACOUSTIC UAUES 


50 


REM 






60 


REM 


THIS PROGRAM 


CALCULATES AND GRAPHS THE NORMALIZED. BACKSCATTERED 


70 


REM 


CROSS SECTION 


FOR ACOUSTICALLY HARD SPHERES OF A SELECTED RANGE 


80 


REM 


OF DIAMETERS D 


ENSONIFIED BY AN INCIDENT PLANE WAUE. 


90 


REM 






100 


REM 






1 10 


REM 






120 


REM 


+ + + + + + + + + + + + + + + + + + + -l- + + + + + -l' + + + + + + + + + + + -E + + -F + + + + + + + + + + + + + -l- + + + -}- + + + -F-l- + + 


130 


! 






140 


I++++ MAIN PROGRAM UARIABLE DECLARATION AND DEFINITION +++++++++++++++++ 


150 






! 


160 






! 


170 


REAL 


Nbs(301 ) 


! THE ARRAY REPRESENTING THE 


180 






! NORMALIZED, BACKSCATTERED CROSS 


190 






! SECTION (NBCS). 


200 






j 


210 


REAL 


K 


! THE WAUE NUMBER, 2*PI*FREQ/C. 


220 






j 


230 


REAL 


A 


! RADIUS OF THE SPHERE. 


240 






I 


250 


REAL 


Ka 


! K*A, THE ARGUMENT OF THE REGULAR 


260 






I AND IRREGULAR SPHERICAL BESSEL 


270 






! FUNCTIONS AND THEIR DERIUATIUES. 


280 






! 


290 


REAL 


Kami n 


! INTEGER OF THE MIN. UALUE OF KA. 


300 






! 


310 


REAL 


K a M a X 


! INTEGER OF THE MAX. UALUE OF KA. 


320 






! 


330 


REAL 


Stp 


! KA INCREMENT. 


340 






! 


350 






! 


360 


•+ + + + + + + + + + + + + + + + + + + + + + -E + + + + + + + + + + + + + + + + + + + + + ++ + + + + + ++ + + + + +++ -I- + + + + + + + + + + 


370 


! 
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380 

390 

400 

410 

420 

430 

440 

450 

4B0 

470 

480 

490 

500 

510 

520 

530 

540 

550 



!++++ MAIN PROGRAM +++++++++++++++++++++++++++++++++++++++++++++++++++++ 



CALL Bcksc tr< Nbs( ♦ ) ,Kamn .Kanax ,9 tp ) ! CALCULATE THE NBS<*). 

j 

CALL Recplt(Nb5( * ) ,Kamn ,Kanax ,5tp ) ! OUTPUT THE NBS(M. 

1 

I +++++++ +++++++++-h++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!++++ MAIN PROGRAM CONCLUSION ++++++++++++++++++++++++++++++++++++++++++ 

I 

PRINTER IS 701 
PRINT CHR$( 12 ) 

END 



1 



+++ + + + +++ + + ++++++ +++ +++++ + +++++++++-}-++ + + + + + +++ + +++++ + 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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5B0 

570 

580 

590 

G00 

B)0 

B20 

B30 

B40 

650 

BB0 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1 030 

1 040 

1050 

1060 

1070 

1080 

1090 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Bcksctr(Nbs( * ) .Kanin .Kanax ,Stp ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE ARRAY NBS( * ) REPRESENTING THE 
REM BACKSCATTERED CROSS SECTION NORMALIZED BY THE CROSS SECTIONAL AREA 
REM OF THE SPHERE. 

REM ++++++ +++++++++++++++++++++++++ +++++++++++++++++++++++++++++++++++ 



!++++ MODULE UARIABLE DECLARATIONS ANO DEFINITIONS +++++++++++++++++++++ 



REAL L 



ORDER OF THE FUNCTION CALCULATED 
IN A SUBROUTINE. IT IS USED 
AS A LOOP INDEX AND TO DESIGNATE 
THE ARRAY ELEMENT FOR THE FUNCTION 
OF CORRESPONDING ORDER. 



REAL Lnax 



MAX. SIGNIFICANT ORDER OF THE SUM 
DETERMINING THE DIFFERENTIAL 
SCATTERING CROSS SECTION (D5CS) 
ANO SCATTERING MODULUS ( SM ) . 



REAL J(51 ) 



ARRAY FOR THE REGULAR SPHERICAL 
BESSEL FUNCTIONS (RSB) OF 
ARGUMENT KA AND ORDER 0 THROUGH 
LMAX. 



REAL 0j(51 ) 



ARRAY FOR THE DERIUATIUES OF THE 
RSB (THE DRSB). 



REAL Y(51 ) 



ARRAY FOR THE IRREGULAR SPHERICAL 
BESSEL FUNCTIONS (ISB) OF 
ARGUMENT KA AND ORDER 0 THROUGH 
LMAX. 



REAL Dy(51 ) 



ARRAY FOR THE DERIUATIUES OF THE 
ISB (THE DISB). 



REAL P(51 ) 



ARRAY FOR THE LEGENDRE 
POLYNOMIALS OF ARGUMENT -1. 



REAL Rsn 



THE REAL COMPONENT OF THE 
SCATTERING MODULUS. 



REAL Isn 



THE IMAGINARY COMPONENT OF THE 
SCATTERING MODULUS. 



REAL Fac 



FACTOR COMMON TO RSM ANO ISM. 



++++ NOTE ++++ 

THE ARRAY DIMENSIONS (30 + 21 ) NEED 
TO BE INCREASED TO ACCOMODATE 
KAMAX>30. 



97 



] ] 00 

1110 

1 1 20 

1 130 

1 140 

1 150 

1 1G0 

1 170 

1 1 80 

1 1 90 

1200 

1210 

1220 

1230 

1 240 

1250 

12B0 

1270 

1280 

1290 

1300 

1310 

1320 

1 330 

1340 

1350 

1 5B0 

1370 

1380 

1390 

1400 

1410 

1420 

1 430 

1 440 

1450 

14B0 

1 470 

1480 

1 490 

1500 

1510 

1520 

1530 

1540 



! 

I 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ CALCULATING THE NBCS ARRAY +++++++++++++++++++++++++++++++++++++++ 

I 

I 



PRINTER IS 1 
INPUT Kanin, Kanax 
Kanin="INT ( Kanin > 

Kanax*I NT( Kanax > 

Sip-. 1 

FOR Ka=Kanin+Stp TO Kanax STEP Stp 
Lnax= INT ( Ka+2 1 > 

CALL Drsb( Lnax ,Ka ,Dj( O ,J( ♦ ) ) 

CALL Di sb( Lnax ,Ka ,Dy( ♦ > ,Y( * ) > 

CALL Legl 80( Lnax ,Ka ,P( ♦ ) ) 

Rsn = 0 
I sn = 0 

FOR L=0 TO Lnax 

Fac = (Dj(L >*(2*L+1 )*P(L)>/(Dj( 
Rsn=R5n+(Fac*Dy(L ) ) 

Isn*=l5n + (Fac*Dj(L) ) 

NEXT L 



! DETERMINE DESIRED RANGE OF KA 
! UALUES. 



! CALCULATE THE DRSB OF ARGUMENT KA 
! AND ORDERS 0 THROUGH LMAX. 

! CALCULATE THE DISB OF ARGUMENT KA 
I AND ORDERS 0 THROUGH LMAX. 

! CALCULATE P(0 OF ARGUMENT ~1 AND 
! ORDERS 0 THROUGH LMAX. 

! CALCULATING THE REAL AND IMAGINARY 
! COMPONENTS OF THE SCATTERING 
! MODULUS. 

L )»Dj(L > + Dy(L >*Dy(L ) ) 



! CALCULATING THE ARRAY REPRESENTING 
! THE NBCS. 



Nbs( (Ka -Kanin >/Stp ) = 4*(R5n*Rsn+l5n*Isn )/(Ka*Ka> 



NEXT Ka 

! 

I 

! ++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++ +++++++++ 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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1550 

15G0 

1570 

1580 

1590 

1G00 

1G10 

1 G20 

1G30 

1G40 

1G50 

1GG0 

1G70 

1G80 

1G90 

1700 

1710 

1720 

1730 

1740 

1750 

17G0 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

18G0 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Drsb ( Lnax ,Ka ,D J ( ♦ ) , J ( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE REGULAR SPHERICAL 
REM BESSEL FUNCTIONS (DRSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX 

REM THE CALCULATION OF THE DERIUATIUE INUOLUES THE RSB’S OF THE SAME 

REM AND PREUIOUS ORDER. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++ CALCULATING THE DERIUATIUES OF THE RSB +++++++++++++++++++++++++++ 



CALL Rsb(Lnax ,Ka ,J( * ) ) 

RAD 

Dj(0) = (-J(0) + COS(Ka) )/Ka 
FOR L=1 TO Lnax 

Dj(D — ((L+1 )/Ka )*J(L ) + J(L-l ) 
NEXT L 



CALCULATE THE RSB UALUES J(0. 

ANGLE UNITS FOR DRSB CALC’NS. 
CALCULATE THE INITIAL DRSB DJ(0). 
CALCULATE THE REMAINING DRSB. 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



SUBEND 



j 



++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++ +++++++ 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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1870 

1880 

1890 

1 900 

1910 

1 920 

1 930 

1940 

1950 

1 960 

1970 

1 980 

1 990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

24 1 0 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Rsb( Lnax ,Ka , J( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE REGULAR SPHERICAL BESSEL FUNCTIONS 
REM (RSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION USES THE CONTINUED FRACTION APPROACH. 

REM + + +++4--I- + + + + + + + + + ++ + + + + + + + + + ++ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + -I- + + 



!++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 

! 

t 



REAL 


R1 


! THE EVOLVING RATIO OF J(0)/J(L). 
1 


REAL 


Nunr 


! THE EUOLUING NUMERATOR FACTOR USED 
! IN CALCULATION OF RL. 

1 


REAL 


Denr 


! THE EVOLVING DENOMINATOR FACTOR 
! USED IN CALCULATION OF RL. 


REAL 


Scexp 


1 

! EU0LUIN6 SCALING EXPONENT FOR UERY 
! LARGE OR UERY SMALL UALUES OF RL. 

1 


REAL 


N f 1 ag 


! 0 FLAG FOR NUMR. 
1 


REAL 


Df lag 


! 0 FLAG FOR DENR. 
1 


REAL 


Bnn 


! EUOLUING TERM USED IN THE 
! CALCULATION OF NUMR AND DENR. 
1 


REAL 


Bine 


! INCREMENT USED IN THE EVOLUTION 
! OF BNM. 

1 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
*♦- + + + CALCULATING THE RSB + + + ++ + + + + + + + + +++ + + + + + + + + + + ++++ + + + + +++ + + + + + + +++ 



RAD 



J(0 ) = ( SIN(Ka ) )/Ka 



ANGLE UNITS FOR CALC'NS. 

CALCULATE THE INITIAL *** 
« * * RSB , J ( 0 ) . ♦ ♦ ♦ 



FOR L=1 TO Lnax 
Rl=l .0 
rJunr = 0 . 

Denr=0 . 
Scexp=0 . 

N f 1 ag=*0 . 

Df lag-0. 

Bnm= 1 . 0/Ka 
B 1 nc-2 . 0/Kd 



CALCULATING THE RATIO *♦* 
J(0 )/J(L). 

INITIALIZING RL. 

INITIALIZING NUMR. 
INITIALIZING DENR. 
INITIALIZING SCEXP. 
INITIALIZING NFLAG. 
INITIALIZING DFLAG. 
INITIALIZING BNM. 

CALCULATE BINC. 
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2420 

2430 

2440 

2450 

24G0 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 



FOR Ordcnt-1 TO L 
Bnm=Bnm+Bi nc 
IF Nflag=0 THEN 
Nunr*Bnm-Numr 
IF Numr<>0 THEN 
Rl^Rl ♦Numr 
Numr= 1 . 0/Numr 

IF ABS(R1 )>10 
ELSE 
Nf lag=-1 
END IF 
ELSE 

Nf iag^O 
END IF 
NEXT Ordcnt 



WHILE NunrODenr OR 
Bnn=Bnm+Binc 

IF Nflag=0 THEN 
Numr«Bnm-Numr 
IF Nunr<>0 THEN 
R1=R 1 *NLimr 
Numr= 1 . 0/Numr 

IF ABS(R1 )>10 
ELSE 
Nf 1 ag=> 1 
END IF 
ELSE 
Nf lag=0 
END IF 

IF Dflag=0 THEN 
Denr“Bnm-^enr 
IF Denr<>0 THEN 
Rl-Rl/Denr 
Denr=1 . 0/Denr 

IF ABS(R1 )>10 
ELSE 

D f 1 ag= 1 
END IF 
ELSE 

Df 1 ag=0 
END IF 

END WHILE 



I * CALCULATE THE L “NAKED" 

I ♦ NUMR FACTORS OF J<0)/J(L). ♦ 

( INCREMENT BNM. 

! CHECK NUMR=0 FLAG NOT SET; 

! UPDATE NUMR IF NUMRO0. 

! UPDATE RL IF UPDATED NUMRO0. 

! PREPARE NUMR FOR NEXT EVOLUTION. 
! CHECK SCALING OF UPDATED RL . 

250 THEN CALL Seal e< R1 .Scexp ) 



! IF NFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



I ♦ CALCULATE THE REMAINING NUMR * 
I ♦ AND DENR TERMS OF J(0)/J(L). ♦ 
! NOTE WHEN NUMR-DENR THEN 
I RL=J(0)/J(L). 

Nflag=1 OR Dflag=1 

! INCREMENT BNM. 

! 

! CHECK NUMR-0 FLAG NOT SET; 

! UPDATE NUMR IF NUMRO0. 

! UPDATE RL IF UPDATED NUMRO0. 

I PREPARE NUMR FOR NEXT EVOLUTION. 
! CHECK SCALING OF UPDATED RL. 

250 THEN CALL Scale( R1 , Scexp ) 



! IF NFLAGO0 THE PREUIOUS STEPS ARE 
I DEFERRED FOR ONE INCREMENT OF BNM. 

! 

! CHECK DENR=0 FLAG NOT SET; 

I UPDATE DENR IF DENRO0. 

I UPDATE RL IF UPDATED DENRO0. 

! PREPARE DENR FOR NEXT EVOLUTION. 

! CHECK SCALING OF UPDATED RL . 

250 THEN CALL Scale( R1 , Scexp ) 



! IF DFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



! 
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2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

31 10 



I 



J(L) = J(0)/(RlM0''(Scexp)) 
NEXT L 



I J(L) DETERMINED. 

I 



! 

! 

+ + + + + + + + + + + + + + + -I-+ + + + +++++ + ++++ + + ++++ + + ++++ + + + + + ++++ ++++++++++ + + + + + +++ + + 

++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



SUBEND 

j 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

! 
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3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Sea 1 e( R1 ,Scexp ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN THE 
REM ONGOING CALCULATION OF RL. 

REM + +++ ++++ ++++ + +++++ + + + + +++++ ++++++++ + + + + + +++ + + + + ++++ + + + + + + + + -I- + + +++ + 

! 

I++++ SCALING RL +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



IF Rl>10''250 THEN 
Rl»Rl*10M-250) 
Scexp=Scexp+250 
END IF 

IF RK10''(~250) THEN 
R1=R1*10"250 
Scexp“Scexp“250 
END IF 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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3430 

3440 

3450 

34G0 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 



j $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Di5b(Lnax,Ka.Dy(* ).Y(0) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE IRREGULAR SPHERICAL 
REM BESSEL FUNCTIONS OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION OF THE DERIUATIUE USES THE ISB’S OF THE SAME AND 
REM PREVIOUS ORDER, 

REM +++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++ +++++++ 

I 

!++++ CALCULATING THE DERIVATIVES OF THE ISB +++++++++++++++++++++++++++ 



CALL Isb(Lnax ,Ka.Y( * ) ) 

RAD 

Dy(0)-(“Y(0) + SIN(Ka) )/Ka 
FOR L-) TO Lnax 

Dy(L)=-((L+l )/Ka )*Y(L )+Y(L-1 ) 
NEXT L 



! CALCULATE THE ISB VALUES. 

! 

I ANGLE UNITS FOR CALC^NS. 

! THE INITIAL ISB DERIVATIVE. 

I THE REMAINING ISB DERIVATIVES, 



j 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



SUBEND 



++++++++++++++ ++++++++++++++++++++++++-h+++++++++++++++++++ +++++++++++++ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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3750 

37G0 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3BG0 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

39G0 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB I sb( Lnax ,Ka , Y( * ) ) 

REM + + + + +4- + + + -H- + + + + + + + + + + + + + + + + + ++ + + + + + ++ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 

REM THIS MODULE CALCULATES THE IRREGULAR SPHERICAL BESSEL FUNCTIONS 

REM OF ARGUMENT KA AND ORDER L THROUGH LMAX . THE CALCULATION IS 
REM ACCOMPLISHED BY FORUARD RECURSION. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



!++++ CALCULATING THE ISB ++++++++++++++++++++++++++++++++++++++++++++++ 



RAD 

Y(0 ) = -COS(Ka )/Ka 

Y( 1 )-( Y(0 )-SIN(Ka ) )/Ka 

FOR L=2 TO Lmax 

Y(L) = ((2*L-1 )/Ka)*Y(L-1 )-Y(L-2) 
NEXT L 



! ANGLE UNITS FOR CALC’NS. 

! INITIAL ISB’S CALCULATED. 

I REMAINING ISB’S CALC’D. 



I 



! 

1 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I + + + + MODULE CONCLUSION + + + + + + + + + + + + + + + + +++ + + + + + + + + + + + + + + + + ■}■ + + + + + + + + + + + + 

I 



SUBEND 



j 



+++++++++++++++++++H-++ +++++++++++++++++++++++++++++++++++++++ ++++++++++ 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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4050 

4060 

4070 

4080 

4090 

4100 

4 110 

4120 

4150 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Leg 1 80( Lnax ,Ka ,P( * ) ) 

REM ++++++++f++++++++++++++++++4+++++ ++++++++++++++++++++ +++++++++++++ 

REM THIS MODULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT -1 AND 
REM ORDER 0 THROUGH LMAX . CALCULATION IS ACCOMPLISHED BY FORWARD 
REM RECURSION. 

REM +f++++++++++++ +++++++++++++++++++++++++++++++++++++++ +++++++++++++ 

j 

!++++ CALCULATING THE LEGENDRE POLYNOMIALS +++++++++++++++++++++++++++H-+ 

! 

X=-l ! ARGUMENT OF THE LEG. POL’S. 

I 

P(0)=1 ! INITIAL LEG. POL’S. CALC’D. 

P( 1 )=X 

FOR L=2 TO Lnax ! REMAINING LEG. POL’S. CALC’D. 

P(L ) = ( 2-( 1 /L ) )*X*P(L-1 )-( l-( 1 /L ) )*P(L-2 ) 

NEXT L 

I 

!+++++++++++++f+++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



SUBEND 



! 



++ +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 






106 



4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4850 

4840 

4850 

4860 

4870 

4880 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Recp 1 1 ( Nbs ( ♦ ) .Kanin .Kanax ,S tp > 

REM + + + + + + + + + + + + -I- + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ++ + + + + + + + + + + + + + + + 

REM THIS MODULE MAPS THE CALCULATED SCATTERING MODULUS ARRAY SM( O ON 

REM A SEMI-LOG PLOT. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

!++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS ++++++4 ++++++++++++++ 

I 

REAL Xgunax ! THE MAXIMUM ABSCISSA UALUE IN 

■ ! GRAPHICS DISPLAY UNITS. 

! 

REAL Ygunax ! THE MAXIMUM ORDINATE UALUE IN 

! GRAPHICS DISPLAY UNITS. 



1 +++++++++++ ++++++T+++++++++++++++++++++++++4-+++++++++++++++++++++++++++ 

!++++ INITIALIZING THE PLOTTER +++++++++++++++++++++++++++++++++++++++++ 
GINIT 

PLOTTER IS CRT .“INTERNAL" 

GRAPHICS ON 
GCLEAR 

j 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
!++++ GRAPH CONSTRUCTION AND LABELLING +++++++++++++++++++++++++++++++++ 



Xgunax-100*MAX( 1 .RATIO ) 

Ygunax=100*MAX( 1 ,1/RATIO) 

! 

CSIZE 4 
LINE TYPE 1 

! 

LORG 4 ! ♦♦ LABELLING THE ABSCISSA AXIS 

MOUE Xgunax/2 ,0 
LABEL "K*A" 

DEG ! ** LABELLING THE ORDINATE AXIS 

LDIR 90 
LORG 6 

MOUE 0 ,Ygunax/2 

LABEL "NORM’D BCKSCTRG CROSS SEC*N" 

LDIR 0 

! ** GRAPH CONSTRUCTION ** 

UIEUPORT 10 .Xgunax-5 . 10 ,Ygunax-5 ! DEFINES THE TOTAL PLOTTING AREA 

! IN GRAPH UNITS. 

l 

UINDOU) Kani n .Kanax ,0 , 1 . 2 
LINE TYPE 4 

GRID Stp*10.. 1 .Kanin.0.61 ,10..2 
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4830 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

51 10 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5130 

5200 

52 10 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 



CLIP OFF 
CSIZE 3 
LINE TYPE 1 

1 

LORG 6 ! LABELLING THE ABSCISSA SCALE ** 

FOR I=Kamin TO Kanax STEP Stp*20 
MOUE I .0 
LABEL I 
NEXT I 

! 

LORG 8 ! ♦* LABELLING THE ORDINATE SCALE 

FOR 1=0 TO 1 STEP .5 
MOVE Kami n , I 
LABEL I 
NEXT I 

I 

CLIP ON 

1 

j 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



I++++ PLOTTING THE CALCULATED NBCS +++++++++++++++++++++++++++++++++++++ 

• ! 

FOR I=Kanin+Stp TO Kamax STEP Stp 
PLOT I ,Nbs( ( I-Kamin )/Stp ) 

NEXT I 



+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++ PRINTING THE PLOT ++++++++++++++++++++++++++++++++++++++++++++++++ 



DUMP GRAPHICS it701 



j 

I 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



I 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



GRAPHICS OFF 
SUBEND 

I 

j + + + + + + + + + + + + + + + + + + + + + + 4+ + + + + + + + 4- + + + + + 4- + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 4' + 
! 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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APPENDIX D 



ACOUSTIC ARRAY EFFICIENCY PROGRAM 



This program was used to analyze data and determine the calibration 



product E^ Ej. The output of this program was also used to construct 



Figures 13, 14, and 15. 



The power supply mentioned in the Gain30 and Gain35 subprograms 



was the HP3314A function generator. All data contained in the data 



subprograms are repeated in Appendix E. Data units are available in 



Appendix E with their respective data 



10 

20 ! 

30 REM + + + 4- + + + -F + + -f- + + + -f + -f-»--^ + -^-l--l--l--F-!- + -«- + -F-l-+-»- + + + + -l- + + + + + + + ++ + 4- + + + + -}- ++ + + + + + + + + + -F 

40 REM ACOUSTIC ARRAY EFFICIENCY 

50 REM, 

G0 REM THIS PROGRAM CALCULATES THE EFFICIENCY AND CENTERLINE GAIN OF 

70 REM AN ACOUSTIC ARRAY BASED ON THE RETURN FROM AN ACOUSTICALLY HARD 

80 REM SPHERE OF DIAMETER D ENSONIFIED BY AN INCIDENT PLANE UIAUE OF 

90 REM FREQUENCY F. 

100 REM 



130 REM 
1 40 I 
150 !++ 



!++++ UARIABLE DECLARATIONS AND DEFINITIONS ++++++++++++++++++++++++++++ 




160 
170 
180 
1 90 
200 
210 
220 
230 
240 
250 



REAL Urnsn 



REAL ’>5P 



MEAN SQUARE UOLTAGE RECEIUED IN 
THE ABSCENCE OF A TARGET AND 
MEASURED AFTER PRE-AMP. 



MEAN SQUARE UOLTA0E RECEIUED f^POM 
THE TARGET SPHERE AND MEASURED 
AFTER PRE-AMPLIFICATION, 
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260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

570 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 



REAL ^nst 
REAL R 



REAL Et 



REAL Go 



MEAN SQUARE UOLTAGE TRANSMITTED. 

THE RANGE FROM THE PLANE OF THE 
ARRAY SPEAKER DIAPHRAGMS TO THE 
CENTER OF THE TARGET SPHERE AND TO 
THE SPECTRUM ANALYZER MICROPHONE. 

THE EFFICIENCY OF CONUERSION 
FROM TRANSMITTED ELECTRICAL POWER 
TO TRANSMITTED ACOUSTICAL POWER. 

THE CENTERLINE GAIN OF THE ACOUSTIC 
ARRAY. 



REAL Etgo 



ET*GO 



REAL At 



THE ATTENUATION COEFFICIENT. 
••ALPHA" . 



REAL C .Cnoi 5 1 



THE SPEED OF SOUND IN HUMID AIR. 



REAL Efreq 

REAL K 
REAL A 
REAL Ka 

REAL Aa 
REAL Gpre 
REAL Nbs 

REAL Atgt 



THE SELECTED FREQUENCY OF THE 
ACOUSTIC ARRAY. 

THE WAVE NUMBER. 2*PI *EFREO/C . 

THE RADIUS OF THE TARGET SPHERE. 

K*A, THE ARGUMENT OF THE REGULAR 
AND IRREGULAR SPHERICAL BESSEL 
FUNCTIONS AND THEIR DERIUATIUES. 

THE AFERATURE AREA OF THE ARRAY. 

THE PRE-AMP’S ELECTRICAL GAIN. 

THE UARIABLE REPRESENTING THE 
NORMALIZED. BACKSCATTERED CROSS 
SECTION (NBCS). 

CROSS SECTIONAL AREA OF THE 
TARGET SPHERE. 



REAL Sasph 



SURFACE AREA OF AN IMAGINARY 
SPHERE CENTERED ON THE ARRAY 
AND HAUING A RADIUS EQUAL TO 
THE RANGE TO THE CENTER OF THE 
TARGET SPHERE. 



REAL Prai 



THE RATIO OF RECEIVED TARGET 
ELECTRICAL POWER TO THE 
TRANSMITTED ELECTRICAL POWER. 



1 10 



790 

800 

810 

820 

830 

840 

850 

850 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1 080 

1070 

1080 

1090 

1 100 

1110 

1 120 

1 130 

I 1 40 

1 150 

1 150 

1 170 

1 180 

1 190 

1200 

1210 

1220 

1230 



REAL £r 



! THE EFFICIENCY OF CONUERSION 
I FROM RECEIVED ACOUSTICAL POWER 
! TO RECEIVED ELECTRICAL POWER. 

! 

} 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

(++++ MAIN PROGRAM +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

I 

CALL Ini t ( Vnsr , Vnsn ,Vns t ,R ,Etgo ,At , A ,Ka , Aa ,Gpre ) 

! INPUT AND CALCULATE REQUIRED 
} PARAMETERS. 

CALL Bcksctr(Ka.Nbs ) ! CALCULATE THE NBCS . 

PRINT "NBS=‘*;Nbs 

! 

Atgt=PI*A*A 
PRINT "ATGT=‘‘ ;Atgt 
Sasph=4 ♦PI *R*R 

Prat = ( Vnsr — Vnsn )/ ( Gpre*Gpre* Vns t ) 

Er=(Prat )♦ ( Sasph^ Sasph/E t go )♦ ( EXP( 2 ♦ At ♦R )/( Aa*Atgt ) )/Nbs 

! 

PRINT “Er^-iEr 
PRINT “EtGo='’ ;Etgo 
PRINT "Go=’‘ jEtgo/Er 
PRINT “EtEr^- ?Er*Er 

) 

! +++ +++ + ++++ + +++ + + + + -+ + + f ++ +++ + + + ++++ +++++ +++ + +++ + + + ++ +++++++++ + + ++++ + + + 

I 

•++++ MAIN PROGRAM CONCLUSION ++++++++++++++++++++++++++++++++++++++++++ 

! 

PRINT CHR$( 12 ) 

PRINTER IS 1 
PRINT 

PRINT "END" ! ADVISES THE USER OF THE MAIN 

PRINT ! PROGRAM’S CONCLUSION. 

END 

! 

}+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +-!' + + + + + + + + ++ + + + + + + + + + + -♦- + + + + +++ + + 






1240 

1250 

12G0 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

13G0 

1370 

1380 

1390 

1400 

1410 

1 420 

1 430 

1 440 

1 450 

I4G0 

1 470 

1 480 

1 490 

1500 

1510 

1520 

1530 

1540 

1550 

15G0 

1570 

1580 

1590 

1G00 

1 G10 

1G20 

1G30 

1G40 

1G50 

1GG0 

1G70 

1G80 

1G90 

1700 

1710 

1720 

1730 

1740 

1750 

I7G0 

1770 

1780 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Init(Un5r,Un5n,Um5t ,R,Etgo ,At .A,Ka ,Aa ,Gpre) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE DETERMINES THE UALUES OF THE PARAMETERS REQUIRED FOR 
REM THE CALCULATION OF Er. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL Ch ! A CONTROL UARIABLE USED FOR 

! CONTROLLING SELECTION OF THE 
I CORRECT GAIN DATA SETS. 

! 

REAL Ur(G0) ! THE ARRAY REPRESENTING THE 

1 MEASURED TARGET RETURN UOLTAGES. 



REAL J ,Jt ,Jn . Jtn 

REAL Tc 
REAL Patm 

REAL Rh 

REAL Tgttn 

REAL Dia 



THE NUMBER OF ARRAY ELEMENTS MINUS 
ONE AND A LOOP CONTROL UARIABLE. 

THE AMBIENT TEMPERATURE IN DEG. C. 

THE ATMOSPHERIC PRESSURE MEASURED 
IN MILLIBARS. 

THE RELATIUE HUMIDITY MEASURED IN 
PERCENT. 

THE MEASURED TIME FOR THE ACOUSTIC 
SIGNAL TO MAKE THE ROUND TRIP 
BETWEEN THE ARRAY AND THE TARGET. 

THE MEASURED DIAMETER OF THE 
TARGET SPHERE. 



REAL Nuspkr 



THE NUMBER OF SPEAKERS IN THE 
ARRAY. 



REAL Diaspkr 



THE AUERAGE DIAMETER OF AN ARRAY 
SPEAKER. 



REAL Ut(G0) 



THE ARRAY REPRESENTING THE 
MEASURED TRANSMISSION UOLTAGES. 



REAL Noise(G0) 



THE ARRAY REPRESENTING THE 
MEASURED NOISE RETURN UOLTAGES. 



REAL Tnoi5e(B0 ) 



REAL Ten 



THE ARRAY REPRESENTING THE 
MEASURED TRANSMISSION UOLTAGES FOR 
NOISE DETERMINATION. 

Tc FOR NOISE DETERMINATION. 
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REAL Rhn 



Rh FOR NOISE DETERMINATION. 



1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1 930 

1 940 

1950 

19G0 

1 970 

1980 

1 990 

2000 

2010 

2020 

2030 

2040 

2050 

20G0 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

21G0 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

22G0 

2270 

2280 

2290 



REAL Patnn 



Patn FOR NOISE DETERMINATION. 



REAL Efreqn 
REAL Cn 
REAL Atn 
REAL Onstn 
REAL Rspkr 



E f reqn=E f req . 

C FOR NOISE DETERMINATION. 

At FOR NOISE DETERMINATION. 

Unst FOR NOISE DETERMINATION. 

THE AVERAGE RADIUS OF AN ARRAY 
SPEAKER. 



I ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++ DETERMINING SYSTEM PARAMETERS +++++++++++++++++++++++++++++++++++++ 

I 

I 

PRINTER IS 1 

PRINT "WHICH TARGET DATA SET 
PRINT "DO YOU WISH ANALYZED?" 

PRINTER IS 701 
INPUT Sphere 

PRINT “TARGET NO. Sphere 
PRINT 

I 

IF Sphere=l THEN 
Ch=3.0 

CALL Gaino(Ch .Etgo ) 

CALL Sphere 1 (Ur ( * ) ,J ,Tc ,Patn ,Rh,Tgttn,Dia ,Efreq .Nuspkr .Diaspkr ,Gpre ) 
CALL Trans50( Jt .Ut ( ♦ ) ) 

CALL Noise30(Jn,Jtn,Noise(* ),Tnoise(* ), Ten, Rhn, Patnn, Efreqn) 

END IF 

! 

IF Sphere=2 THEN 
Ch=3.5 

CALL Ga i no( Ch ,Et go ) 

CALL Sphere2(Ur( ♦ ) ,J ,Tc ,Patn ,Rh ,Tgttn ,Dia ,Efreq , Nuspkr .Diaspkr ,Gpre ) 
CALL Trans35( Jt ,Ut( M ) 

CALL Noise 35 ( Jn ,Jtn ,Nolse ( * ) ,Tnoise( * ) ,Tcn ,Rhn , Patnn .Efreqn ) 

END IF 

I 

IF Sphere-3 THEN 
Ch=3.5 

CALL Gai no( Ch ,Et go ) 

CALL Sphere3(Ur( • ) ,J ,Tc ,Patn ,Rh .Tgttn ,Dia .Efreq .Nuspkr .Diaspkr ,Gpre ) 
CALL Trans55( J t ,Ut ( * ) ) 

CALL Noise35(Jn,Jtn,Noise(» ),Tnoise( * ), Ten, Rhn, Patnn, Efreqn) 

END IF 
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! 



2300 

2310 

2320 

2330 

2340 

23B0 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 



IF Sphere*4 THEN 
Ch=3.5 

CALL Gaino( Ch .Etgo ) 

CALL Sphere4( Ur( ♦ ) ,J ,Tc ,Patn ,Rh,Tgttn ,Dia ,Efreq ,Nuspkr .Diaspkr ,Gpre ) 
CALL Trans35( J t ,Ut ( ♦ ) ) 

CALL Noi5e35( Jn,Jtn,Noise( ♦ ) ,Tnoise ( * ) ,Tcn,Rhn ,Patnn ,Efreqn) 

END IF 



CALL Uol t sns ( Ur ( ♦ ) , J , Vnsr ) 

PRINT "UNSR^" ;Unsr 

CALL Atten(Tcn,Rhn,Patnn ,Cn ,Atn .Efreqn) 

CALL Uolt5ms(Tnoi5e( ♦ ) ,Jtn ,Unstn ) 

CALL Uol t sns( Noi se( ♦ ) , Jn ,Unsn ) 

CALL Uolt5ns(Ut( ♦ ) . Jt ,Unst ) 

PRINT "UMST=" ;Unst 

I 

A=Dia/2 
PRINT ’•A='^A 

CALL Range(Tgttn ,A ,Tc ,Rh,Patn ,C ,At .Efreq ,R) 

PRINT “C=“;C 

PRINT "AT=”;At 

PRINT “ATTEND" ;EXP(At*R ) 

PRINT ”R=";R 

Unsn=Umsn*EXP( 2*R* ( At n-At ))*Un5t/Unstn 

PRINT "UMSN=" ;Unsn 

Ka=2*PI*Efreq*A/C 

PRINT "KA-“;Ka 

Rsp kr=Di aspkr /2 

Aa=Nuspkr*PI*Rspkr*Rspkr 

PRINT ’‘AA^’^Aa 



I-+ + + + + + + + + H 



!++++ NODULE CONCLUSION ++++■ 



K -}-+++++++++++++++ 



K+++++ ++++++++++ 



! 

SUBEND 



! 

! +++ + + + + +++++-I--I- ++ +++ + +++++++++ + + + ++++ + + +++++++ ++++++++ ++ + + ++++ 
I 

! $$$$$$$$$$ 5 $$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$ 
I 
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2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

28!0 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

31 10 

3120 

3 1 30 

3140 

3150 



SUB 6ai no( Ch ,Etgo ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE DETERMINES THE UALUE OF EtGo. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++ 
I 

I++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL Urns 



ROOT MEAN SQUARE UOLTAGE MEASURED 
IN GAIN DETERMINATION. 



REAL Ir 



MEASURED INTENSITY AT R; 
Ptgt /At g t . 



REAL Za 



ACOUSTIC ARRAY IMPEDANCE. 



^++++++++H 



+ + + + + + + + + + + + + + + + + + + + + + + + + -r + + ++ + + + + + + + + + + + + + + + + 



++++ DETERMINATION OF EtGo + 



++++++++T++++++ +++++++++++++ 



I 



IF Ch=3.0 THEN CALL Gai n30( Urns ,Tc ,Rh ,R , Pat n , Ir ,E freq , Za ) 

IF Ch=3.5 THEN CALL Gain35( Urns ,Tc ,Rh ,R ,Patn , Ir ,Efreq ,Za ) 

! OBTAINING UALUES FROM APPROPRIATE 
! DATA. 

! 

CALL At ten( Tc ,Rh ,Patn ,C , At .Efreq ) ! CALCULATING AT. 

Pt=Urns*Urns/Za 

Etgo=Ir*4*PI*R*R»EXP( At»R )/Pt I CALCULATING EtGo. 

I 



1 



+++++H 



I- + + + 



+ + + + + + + + + + + ++ + + + + + + + - 1 -+ + + + + ^“ -!- + + ->- + + + + + + + + + + + + + ++ + + + 
++++ MODULE CONCLUSION +++++++++++++++++++++++++++++++++++++H 
SUBEND 

++++++++++++++++++++++++++++++++++ ++++++++++++++++T+i‘+ ++++++ +++++++++++ 

$ss$s$ss$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$ 



I 



I 
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3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 



I 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Uolt5ns(Uolt( * ) .J ,Uns) 

REM +++++++++++++++++++++++++++++ ++++-K+ +++++++++++++++++++++++++++++++ 
REM THIS MODULE DETERMINES THE MEAN SQUARE UALUE OF THE INPUT VOLTAGE 
REM ARRAY. 

REM ++++++++++++++++++4*+++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++ MODULE VARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 

I 



!REAL VOLT( * ) 

REAL Volts 
REAL Vavg 
REAL V 

I REAL Vns 



! 

1 THE ARRAY REPRESENTING THE INPUT 
VOLTAGE ARRAY. 

THE SUM OF V0LT(O’S ELEMENTS. 

THE AVERAGE VALUE OF V0LT(O. 

THE DIFFERENCE BETWEEN VOLTd > 

AND VAVG. 

1 THE MEAN SQUARE VALUE OF VOLT(*>. 






h++ +++++++++++++++++ +++++++++++ 



! ++++ DETERMINATION OF VMS ++T++++++++-*-+-rT-r4-++^-^++++++++ ++++++++++++++-i' + 

I 

I 

Voits=0 
FOR 1=0 TO J 

Volt5=Volts+Volt(I ) 

NEXT I 

Vavg=Volts/( J+1 ) 

FOR 1=0 TO J 
V=Volt( I )-Vavg 
Vns=( V*V)+Vns 
NEXT I 

Vns=Vns/ ( J+ 1 ) 

! 

! 

}+ + + + + ++ + + + + +-i- + ++ + '*- + -*--^4.-;- + + + 4. + + 4-+ + 4. + + + + + 4. + -(.-r + + +-r + + + + + + + 'r-'- + i' + + + + + + + i--r-^ + + + + 
I 

I + + + + MODULE CONCLUSION + + + + + + + + ^ + + + •** ++4“ + + -^ + + + + -r + + + *^-''++ +++ + 'f + + + -+' + + + ++++ + 

I 

SUBEND 

I 

!+ + + + + + + + + i- + + + + + + -t-+ + + + + + -f- + + + -i- + + + + + + + + r + + -^-^ + + -K + -t‘->' + -K+-J- + + + + + + + r-^ + + + ^-rr + -!- + + -l' 
I 

I $$$$$$s$$$$$$$s$$ $$$$$$$$$$$$$$$$$$$ $$SSS$S$$SS5$SS$$$S$SS$$SS$SS$SSSSS 
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3B70 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 



I 



I $$$$$$$$$$$$$$ 

SUB Range ( Tgttn,A,Tc , Rh, Pa tn,C, At ,Efreq,R) 

REM THIS MODULE DETERMINES THE AVERAGE PATH ATTENUATION AND THE RANGE 
REM BETWEEN THE ARRAY AND THE TARGET CENTER. 



!++++ DETERMINATION OF R AND AT 



K+++++++ 



Tgt tn-Tgttn/1000 

CALL Atten(Tc,Rh,Patn,C,At ,Efreq) 
R=Tgttn*C/2+A 



AT AND C DETERMINED. 
R DETERMINED. 



++++ MODULE CONCLUSION 
SUBEND 



h++++++++++H 



h^++++++++++ +++++++++ 



h +++++++++++++++ ++++++ 



$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$ 



I 



5930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4 100 

4110 

4120 

4130 

4140 

4 150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 



! 

I $$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$ 
SUB Atten( Tc ,Rh ,Patn .Cnoist ,At .Efreq ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE ATTENUATION USING RELATIVE HUMIDITY 
REM IN PERCENT. ATMOSPHERIC PRESSURE IN MILLIBARS. AND TEMPERATURE 
REM IN DEGREES CENTIGRADE. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



I++++ MODULE VARIABLE DECLARATIONS 

REAL Tk 
REAL Es 
REAL Pratio 

REAL Cdry 
REAL H 
REAL Past 
REAL Tast 
REAL Fn 
REAL Anax 
REAL Frat 
REAL F2 
REAL Acl .Anol 



AND DEFINITIONS +++++++++++++++++++++ 

j 

I 

I THE AMBIENT TEMPERATURE IN DEG. K. 

! 

! THE SATURATION VAPOR PRESSURE. 

I 

! THE RATIO OF WATER PRESSURE IN 
! MILLIBARS TO ATMOSPHERIC PRESSURE 
! IN MILLIBARS. 

I 

! THE SPEED OF SOUND IN DRY AIR. 

I 

I 

! 

! P ASTERISK. 

! 

! T ASTERISK. 

j 

i MAX. FREQ. 

I 

! MAX EXPECTED ATTENUATION. 

I 

! EFREQ/FM 

i 

I FRAT*FRAT 
I 

i ATTENUATION EXPONENT COMPONENTS. 



I + + ++++ + ++++ + + + +++ + + +++ + + + -». + + + + + + + + -I 



h++ + + +-l- + + + + + + + H 



I- + + + + 
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4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 



! ++++ CALCULATING THE ATTENUATION ++++++++-M-+++++++++++++++H 



^ + + + + - 1 - + + 



Tk=Tc+273. 15 



Rh=Rh/100 

Es=10"( 9,4-2353/Tk ) 
Pratio=E5*Rh/(Patrn-Es*( 1+Rh)) 



I 

I 

! CONUERSION OF TEMP IN CELSIUS TO 
! TEMP IN KELUIN. 

! 

! CONUERSION OF RH IN PCT . TO PRATI0. 
1 



Cdry=20.05*SQR(Tk ) 

Cnoi si=Cdry*( 1 + . 14*Pratio) 



I CALCULATION OF THE SPEED OF SOUND. 



I 



H=Pratio*100 

Pas t-Patn/ 1014 

Tast = ( 1 .8»Tc + 492 )/519 

Fn=( 10 + 6600*H+44400*H*H)*Past/(Tast'' .8 ) 

! 

Anax-.0078*Fn*(Tast"( -2.5 ) )*EXP(7.77»( 1-1/Tast ) ) 

! 

Frat =Ef req/Fn 
F2 = Frat *Frai 

I 

Anol=Anax*SQR( .18*. 1 8*F2 + ( 2 *F2/ ( 1+F2 ) )"2 )/304.8 
Acl = l .74*10"(-10 )*Efreq*Efreq 
At = ( Amol+Acl )*LOG( 10 )/ 10 



h + + + + + + + + + +-l- + + + + + + + + + + + -f- + + + + + + + + + + H 



f--^ + + + + + + H 



H+++++++ 



++++ MODULE CONCLUSION + + + + ++++ + ++++++ +++ + + + +++ + +++ ++-?'++ +++++ +++ + + + +++ + 

I 

SUBEND 

j 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

2 $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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4720 

4730 

4740 

47B0 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

51 10 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 



I 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Bcksc tr ( Ka ,Nb5 ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE BACKSCATTERED CROSS SECTION NORMALIZED 
REM BY THE CROSS SECTIONAL AREA OF THE SPHERE. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

I++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL L 



ORDER OF THE FUNCTION CALCULATED 
IN A SUBROUTINE. IT IS USED 
AS A LOOP INDEX AND TO DESIGNATE 
THE ARRAY ELEMENT FOR THE FUNCTION 
OF CORRESPONDING ORDER. 



REAL Lnax 



MAX, SIGNIFICANT ORDER OF THE SUM 
DETERMINING THE NORMALIZED 
BACKSCATTERING CROSS SECTION 
(NBCS ). 



REAL J(51 ) 



ARRAY FOR THE REGULAR SPHERICAL 
BESSEL FUNCTIONS ( RSB ) OF 
ARGUMENT KA AND ORDER 0 THROUGH 
LMAX. 



REAL Dj(51 ) 



ARRAY FOR THE DERIUATIUES OF THE 
RSB (THE DRSB). 



REAL Y(51 ) 



! ARRAY FOR THE IRREGULAR SPHERICAL 
! BESSEL FUNCTIONS (ISB) OF 
I ARGUMENT KA AND ORDER 0 THROUGH 
I LMAX. 



REAL Dy(5l ) 



ARRAY FOR THE DERIUATIUES OF THE 
ISB (THE DISB). 



REAL P(51 ) 

REAL Fac 
REAL Rsn 

REAL Isn 



ARRAY FOR THE LEGENDRE 
POLYNOMIALS OF ARGUMENT COS(180). 

FACTOR COMMON TO RSM AND ISM. 

THE REAL COMPONENT OF THE 
SCATTERING MODULUS, SM. 

THE IMAGINARY COMPONENT OF THE 
SM. 
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5210 

5220 

5230 

5240 

5250 

52G0 

5270 

5280 

5290 

5300 

5310 

5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 

5620 



+ + *♦-+ NOTE + + + + 

THE ARRAY DINENSIONS (30 + 21 ) NEED 
TO BE INCREASED TO ACCOMODATE 
KA>30. 



I ++++++++++++++++++ ++++++++++++++ 
! 

I++++ CALCULATING THE NBCS ++++++ 
Lnax=INT ( Ka+2 1 ) 

'CALL Drsb ( Lmax ,Ka ,D j ( ♦ ) , J ( ♦ ) ) 



CALL Dlsb(Lnax ,Ka ,Dy( ♦ ) ,Y( ♦ ) ) 



CALL Leg180(Lnax ,Ka ,P( * ) ) 

Rsn=0 
I sn = 0 

FOR L=0 TO Lnax 

Fac = (Dj(L)*(2*L+1 )*P( L ) )/ ( D j ( L ) * 
Rsn=^Rsn + ( Fac*Dy ( L ) ) 
Isn=Ism+(Fac*Dj(L ) ) 

NEXT L 

Nb5 = 4*(R5n*Rsrn + Isn*l5n )/(Ka*Ka ) 



++++++++++++++++++++++++++++++++++++ 

++++++++++++++++++++++++++++++++++++ 



! CALCULATE THE DERIUATIUES OF THE 
! REGULAR SPHERICAL BESSEL FUNCTIONS 
! OF ORDER 0 THROUGH LMAX . 

! CALCULATE THE DERIUATIUES OF THE 
! IRREGULAR SPHERICAL BESSEL FUNC’NS 
! OF ORDER 0 THROUGH LMAX . 

! 

! CALCULATE THE ARRAY Pi*), 

! 

! CALCULATING THE REAL AND IMAGINARY 
! COMPONENTS OF THE SM. 

j(L)+Dy(L)*Dy(L)) 



CALCULATING THE NBCS. 



+ + + +++ + + +++ + ++ +++ + + + + + + + + + + + + + + ++++-^ + + + + + + + + +++ +++++ ++++ + + + + + + +++ + + + + + + 



++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



SUBEND 

! 

!+ + + + + + + + + + + + + +. + + + + + + + + + + + + + + + + + + + + + + + + +r + + + + + ^ + + +^ + ^ + + + + -f-+ + + + + + ++ + + + + + + 



$SS$$S$SSSSSSS$SSS$$S$S$$$$S$$$$$SS$$S$S$$SSSS$S5SS$SSS$SSSSSSS$S$S$SSS 
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5B30 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$ $$$$$$$$$ 

SUB Drsb(Lnax .Ka ,Dj( * ) ,J( * ) ) 

REM ++++++++++++++++++++++++++++++++++++++ +++++++++++++1-+ +++++++ ++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE REGULAR SPHERICAL 
REM BESSEL FUNCTIONS (DRSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . 

REM THE CALCULATION OF THE DERIUATIUE INUOLUES THE RSB’S OF THE SAME 

REM AND PREUIOUS ORDER, 

I 

!+ +++ CALCULATING THE DERIUATIUES OF THE RSB +-^ +++-^ ++++++ +++++++++++ ++++ 

} 



CALL Rsb(Lnax ,Ka , J< ♦ ) ) 

RAD 

Dj(0) = (-J(0)+COS(Ka))/Ka 
FOR L=l TO Lnax 

Dj(L)=-((L + l )/Ka )*J(L) + J(L-l ) 
NEXT L 



CALCULATE THE RSB UALUES J(*). 

ANGLE UNITS FOR DRSB CALC’NS. 
CALCULATE THE INITIAL DRSB DJ(0). 
CALCULATE THE REMAINING DRSB. 



++^+ MODULE CONCLUSION ++++++++H 
SUBEND 



h +++++++++++++++++++++ +++++++++ 
hr + + -K + + + + + + + + + + + 'r + + + + + + + + 



+ + + + + + + + + + -•“ + + + + + + + + + + + + + ++ + + + + + + + ++ + + + + + + + ++ + + H 



I- + + + + +H 



$$$$ $$$$$$$$ 2 $ $$$$$$$$$$$$$$$$$$$$$$$$ $$$$ 5 $$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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5950 

59G0 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

61 10 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 

6560 

6370 

6380 



! 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Rsb ( Lnax .Ka , J( ♦ ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE REGULAR SPHERICAL BESSEL FUNCTIONS 
REM <RSB) OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION USES THE CONTINUED FRACTION APPROACH. 

REM ++++++++++++++++++++++++++++++++++++++++ +++++++1-++++++++ ++++++++++ 

! 

!++++ MODULE UARIABLE DECLARATIONS AND DEFINITIONS +++++++++++++++++++++ 



REAL R1 
REAL Nunr 

REAL Denr 

REAL Scexp 

REAL Nflag 
REAL Dflag 
REAL Bnn 



THE EUOLUING RATIO OF J(0)/J(L). 

THE EUOLUING NUMERATOR FACTOR USED 
IN CALCULATION OF RL. 

THE EUOLUING DENOMINATOR FACTOR 
USED IN CALCULATION OF RL. 

EUOLUING SCALING EXPONENT FOR UERY 
LARGE OR UERY SMALL UALUES OF RL. 

0 FLAG FOR NUMR. 

0 FLAG FOR DENR. 

EUOLUING TERM USED IN THE 
CALCULATION OF NUMR AND DENR. 



REAL Bine 



INCREMENT USED IN THE EVOLUTION 
OF BNM. 



! ++++++++++++++H 



I 



h+++++++++rH 



^+++++++++++ ++++++++ 



!++++ CALCULATING THE RSB +++++++++++++++++++++++++++++++H 

I 



“+++++++ 



RAD 

J(0) = (SIN(Ka))/Ka 



ANGLE UNITS FOR CALC’NS. 

♦♦♦ CALCULATE THE INITIAL ♦♦♦ 
♦♦♦ RSB, J(0). 
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G390 

G400 

G4I0 

G4Z0 

G450 

G440 

G450 

G4G0 

G470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 



FOR L=1 TO Lnax 
R1=I .0 
Nunr=0 . 

Denr=0 . 
Scexp=0 . 

N f 1 ag = 0 . 

Of lag=0. 

Bnn=1 .0/Ka 
B i nc=2 . 0/Ka 



CALCULATING THE RATIO *** 
***J(0)/J(L). 

INITIALIZING RL. 

INITIALIZING NUNR. 
INITIALIZING DENR. 
INITIALIZING 5CEXP. 
INITIALIZING NFLAG. 
INITIALIZING DFLAG. 
INITIALIZING BNM. 

CALCULATE BINC. 



FOR Ordcnt=1 TO L 
Bnrn = Bnn + Bi nc 
IF Nflag=0 THEN 
Nunr*Bnrn-Nunr 
IF Nunr<>0 THEN 
R 1 =R 1 *Nunr 
Numr= 1 . 0/Nunr 

IF AB5(R1 )>10"250 THEN CALL 
ELSE 
Nf lag=1 
END IF 
ELSE 

Nf 1 ag=0 
END IF 
NEXT Ordcnt 



WHILE NunrODenr OR Nflag=1 OR Df 
Bnrn=Bnn + Bi nc 

IF Nflag=0 THEN 
Nunr=Bnn-Nunr 
IF Nunr<>0 THEN 
R1 ="R1 *Nunr 
Nunr= I . 0/Numr 

IF ABS(R1 )>10''250 THEN CALL 
ELSE 

N f 1 ag= 1 
END IF 
ELSE 

N f 1 ag = 0 
END IF 

IF Dflag=0 THEN 
Denr=Bnn-Denr 
IF Denr<>0 THEN 
R1=R1/Denr 
Denr= 1 . 0/Denr 



I * CALCULATE THE L “NAKED” 

! ♦ NUNR FACTORS OF J(0)/J(L). * 

! INCREMENT BNM. 

! CHECK NUMR=0 FLAG NOT SET; 

I UPDATE NUMR IF NUMRO0. 

f UPDATE RL IF UPDATED NUMRO0. 

I PREPARE NUMR FOR NEXT EUOLUTION. 
! CHECK SCALING OF UPDATED RL . 

Seal e( R 1 , Seexp ) 



! IF NFLAGO0 THE PREVIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



! * CALCULATE THE REMAINING NUMR * 
I « AND DENR TERMS OF J(0)/J(L). * 
! NOTE WHEN NUMR=DENR THEN 
! RL=J(0 )/J(L ). 
ag=1 

! INCREMENT BNM. 

I 

! CHECK NUMR=0 FLAG NOT SET; 

! UPDATE NUMR IF NUMRO0. 

I UPDATE RL IF UPDATED NUMRO0. 

! PREPARE NUMR FOR NEXT EUOLUTION. 
I CHECK SCALING OF UPDATED RL . 

Seal e( R1 , Seexp ) 



! IF NFLAGO0 THE PREUIOUS STEPS ARE 
! DEFERRED FOR ONE INCREMENT OF BNM. 



I CHECK DENR=0 FLAG NOT SET; 

I UPDATE DENR IF DENRO0. 

! UPDATE RL IF UPDATED DENROO. 

I PREPARE DENR FOR NEXT EUOLUTION. 



124 



6930 

6940 

6950 

6960 

6970 

6980 

6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

71 10 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 



! CHECK SCALING OF UPDATED RL. 

IF ABS<R1 )>10"250 THEN CALL Scdle( R1 .Scexp ) 

ELSE 
Df lag= 1 
END IF 

ELSE ! IF DFLAGO0 THE PREVIOUS STEPS ARE 

Dflag=0 ! DEFERRED FOR ONE INCREMENT OF BNM. 

END IF 

! 

END WHILE 



J(L)=J<0)/<Rl*10"(Scexp)) 
NEXT L 



J(L) DETERMINED. 



!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

j 

SUBEND 

j 

I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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7200 

7210 

7220 

7230 

7240 

7250 

72G0 

7270 

7280 

7290 

7300 

7310 

7320 

7530 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7450 

7440 

7450 

7460 

7470 

7480 

7490 

7500 



I 



I $$$$$$$$$$$ $$$$$$s$$$$$$$$$$$s$$$$s$$ss$$$$$s$$$$$$ss$s $$$$$$ $s$$ss$$ss 

SUB Seal e( R1 ,Scexp ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE SCALES RL TO MAINTAIN NUMERICAL ACCURACY IN THE 
REM ONGOING CALCULATION OF RL . 

I 

I++++ SCALING RL +++++++++++++++++++++++++++++++++++++++++++++++^ 



IF Rl>10''250 THEN 
Rl=Rl*10"(-250) 

Scexp=Scexp+250 
END IF 

IF Rl<]0"(-250) THEN 
Rl=Rin0"250 
Sc ex p= Sc exp ”2 50 
END IF 

j 

j 

I + -^ + + ->- + + + + + + + + + + + + + + -l- + + + + + + '^ + -l- + + + + + + + + + + + + + + + + + + + + + + + + + + + + •!“ + + + + + + + + + + + + + 



++++ MODULE CONCLUSION 
SUBEND 



I- + 4- + + + + + + + H 



I- + + + + + + + + + + +-^ + + ++ + + + T + + + + + + + -f-+ + +++ + TH 



I 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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7510 

7520 

7530 

7540 

7550 

75B0 

7570 

7580 

7590 

7600 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Disb(Lnax ,Ka ,Dy( M ,Y( O) 

REM +++++++++++++++++++++++++++++++++++++++++++++ ++++++ ++++++ 
REM THIS MODULE CALCULATES THE DERIUATIUES OF THE IRREGULAR SPHERICAL 
REM BESSEL FUNCTIONS OF ARGUMENT KA AND ORDERS 0 THROUGH LMAX . THE 
REM CALCULATION OF THE DERIUATIUE USES THE ISB’S OF THE SAME AND 
REM PREUIOUS ORDER. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

f 

•++++ CALCULATING THE DERIUATIUES OF THE ISB +++++++++++++++++++++++++++ 



CALL IsbCLnax ,Ka ,Y( ♦ ) ) 

RAD 

Dy(0 )=(-Y(0 ) + SIN(Ka ) )/Ka 
FOR L=1 TO Lnax 

Dy(L)=-((L+1 )/Ka )»Y(L) + Y(L-1 ) 
NEXT L 



! CALCULATE THE ISB UALUES. 

! ANGLE UNITS FOR CALC’NS. 

! THE INITIAL ISB DERIUATIVE. 

! THE REMAINING ISB DERIUATIUES. 

I 

I 



I +++++++H 



++++++++++++++++++++++++++++++++++++++++++++++++ 



I 

>++++ MODULE CONCLUSION ++++++++++++++++++++++++++++T+++++++++++++++++++ 

I 

SUBEND 

I 



++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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7830 

7840 

7850 

7880 

7870 

7880 

7890 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

7970 

7980 

7990 

8000 

8010 

8020 

8030 

8040 

8050 

8060 

8070 

8080 

8090 

8100 

81 10 

8120 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $s$$$$ $$$$$$$ 
SUB Isbdnax .Ka ,Y( ♦ ) ) 

REM +++++++++++++++++++++++++++++++++++ +++J-+ +++++++++++++ +++++++++++++ 

REM THIS MODULE CALCULATES THE IRREGULAR SPHERICAL BESSEL FUNCTIONS 
REM OF ARGUMENT KA AND ORDER L THROUGH LMAX . THE CALCULATION IS 
REM ACCOMPLISHED BY FORWARD RECURSION. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

•++++ CALCULATING THE ISB ++++++++++++++++++++++++++++++++++++++++++++++ 



RAD ! ANGLE UNITS FOR CALC’NS. 

Y(0)=-COS(Ka)/Ka ! INITIAL ISB’S CALCULATED. 

Y( 1 )-( Y( 0)“SIN(Ka ) )/Ka 

FOR L=2 TO Lnax ! REMAINING ISB’S CALC’D. 

Y(L) = ((2*L~1 )/Ka)*Y(L-l )-Y(L-2) 

NEXT L 

I 

I 

1 +++++++++++++++++++++++++++++++++++++++++++++-^+^+++++++++++++++++++++++ 



!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

I 

!+++ + + + + + + + + + + + + + + + + + + + -l- + + + -|- + + + + + + + + + -^ + + + 'r + + -f- + T + + + -l‘ + + + + + + + + + + + ++ + + + +'r + + + 

i $$$$$$$$$$$$$$S$$$$$$S$$$$$$$$S$$$$$S$SS$$$$S$$$S$$$$$$$$S$SS$$$SS$$S$$ 



128 



8130 

8140 

8150 

8160 

8170 

8180 

8190 

8200 

8210 

8220 

8230 

8240 

8250 

8260 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8360 

8370 

8380 

8390 

8400 

8410 



I $$$$$$$S$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$$$$$$$$$$S£S$$$$SS$$$S$$$S$$$S$S 
SUB Legl80(Lnax .Ka ,P( ♦ ) ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CALCULATES THE LEGENDRE POLYNOMIALS OF ARGUMENT -1 AND 
REM ORDER 0 THROUGH LMAX . CALCULATION IS ACCOMPLISHED BY FORWARD 
REM RECURSION, 

REM ++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++ 

! 

!++++ CALCULATING THE LEGENDRE POLYNOMIALS +++++++++++++++++++++++++++++ 

! 

X=-l ! ARGUMENT OF THE LEG. POL’S. 

j 

P(0)=1 ! INITIAL LEG. POL’S. CALC’D. 

P( 1 ) = X 

FOR L=2 TO Lnax ! REMAINING LEG. POL’S. CALC’D. 

P(L ) = (2“( 1 /L ) )*X*P(L-1 )-( l-( 1/L ) )*P(L“2 ) 

NEXT L 

! 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



++++ MODULE CONCLUSION 
SUBEND 



r +++++++++++++++++++++++++++++ 



h +++++++++++++++++++++++++++++++++++++++++ +++++++++ 



1 



t 
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8420 I 

8430 I $$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$5$$$$$$$ $$$$$$$ 
8440 SUB Gain30( Urns ,Tc ,Rh ,R ,Patn , Ir .Efreq ,Za ) 

8450 REM +++ +++ + + ++++ + + + + + + + + + + + + + ++'♦' +++ + + + + + + +++ +++ + + + + + + + +++++++++++ + +++ + 

8460 REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF EtGo UITH 

8470 REM 3.0U INPUT FROM THE POWER SUPPLY. 

8480 REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8490 I 

8500 !++++ LOADING DATA ++++++++++++++++++++++T++++++++++++++++++++++++++++++ 

8510 ! 

8520 ! 

8530 RESTORE 630 

8540 READ Urns ,Tc ,Rh ,R ,Patn , Ir ,Efreq ,Za 

8550 G30: DATA .426.21.0,50.0,5.30,1021.5,2.9515-3.5000,15.8 

8560 ! 

8570 ! 

8580 ! ++++++++ ++++++++ ++++++++++++++++++++++++++++++++++++++ +++++++++++++++++ 

8590 1 

8600 !++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

8610 I 

8620 SUBEND 

8630 ! 

8640 ! + + + + + +++++++ ++ + -^ + + ++ + + + +++++ + +++++++ + + ++++'♦'++++ + + + ++++++++++ ++ + + + + + +++ + 

8650 I 

8660 « $$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ 
8670 ! 
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8G80 ' 

86 90 I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
8700 SUB Gain35( Urns ,Tc ,Rh ,R ,Patn , Ir ,£freq , 2a ) 

8710 REM ++++++ + + + +++++ +++++ + ++++++ + + +++++++-f + +++ +++ + ++++ ++ ++++++ +++++++ +++ 

8720 REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF EtGo WITH 

8730 REM 3.5U INPUT FROM THE POWER SUPPLY. 

8740 REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8750 « 

8760 ! ++++ LOADING DATA +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8770 ! 

8780 ! 

8790 RESTORE G55 

8800 READ Urns ,Tc ,Rh ,R ,Patn , Ir ,E freq ,Za 

8810 G35: DATA . 523 ,2 1 . 0 .50 . 0 ,5 . 30 , 1 02 1 . 5 .4 . 467E-3 ,5000 . 1 5 . 8 

8820 I 

8850 ! 

8840 I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
8850 ! 

8860 ! ++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

8870 ! 

8880 SUBEND 

8890 ! 

8900 ! +++++++++++ ++++++++++++++++++++++++++++++++++++++T+ ++++++++++++++++++++ 

8910 ! 

8920 ' $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

8930 ! 



8940 
8950 
89B0 
8970 
8980 
8990 
9000 
90)0 
9020 
9050 
9040 
9050 
9080 
9070 
9080 
9090 
9100 
91 10 
9120 
9130 
91 40 
9150 
9180 
9170 
9180 
9190 
9200 
9210 
9220 
9230 



( 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Sphere 1(Dt< * ) ,J ,Tc ,Patn,Rh,Tgttn,Dia,Efreq .Nuspkr.Diaspkr.Gpre ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Er WITH 
REM TARGET NO. 1 EMPLOYED. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

■+++++++++++++++++ +++++++++++++H 
I 
I 



»++++ LOADING DATA 



I-+++ +++++++ 



J = 39 

RESTORE SI 
FOR 1=0 TO J 
READ Dt( I ) 

NEXT I 

READ Tc, Pa in, Rh,Tgttn,Dia,Efreq.Nuspkr,Diaspkr,Gpre 
SI ; DATA 9 . 54 , 9 . 4 . 8 . 32 .8 . 42 ,3 . 88 . . 98 .-2 . 07 ,-4 . 9 ,-7 . 24 ,-8 . 87 

DATA -9.84,-9.45,-8.32,-8.41 .-3 . 83 .- . 89 .2 . 1 4 .4 . 94 .7 .28 .8 . 88 
DATA 9. 8. 9. 38, 8. 22. 8. 28. 3. 87,. 73, -2. 5. -5.1 1 ,-7.4, -8. 98 
DATA -9.84 ,-9.39 ,-8.2 ,-8.21 ,-3.81 ,-.85,2.38.5.12,7.38,8.91 
DATA 20.2,1015.8.52.8.29.91 ,. 2545 .5000 , 1 9 ,. 0752 , 1 1094 

i 

I 

! +-»► + + + + + + + + + + + + + ++ + + + + + + + + + + + + + + + + + + + + + + +-»■ + + + + -^ + + + + + + + + -^ + + + + + + + + + + + + + + + + 
I 

!++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 



I 



SUBEND 



9240 ! + + + ++++ + +++ + +++ + ++++ + +++ + ++++++++++ +++++ ++++++++-^++-^ ++++++ + ++++ +++ +++ + + 

9250 ! 

9250 • $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$s$$$$$$$S $$$$$$$$$$$$$$$$$$$$ $$$$$$$$ 

9270 « 
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9280 

9290 

9300 

9310 

9320 

9330 

9340 

9350 

9360 

9370 

9380 

9390 

9400 

9410 

9420 

9430 

9440 

9450 S 

9460 

9470 

9480 

9490 

9500 

9510 

9520 

9530 

9540 

9550 

9560 

9570 

9580 

9590 

9600 

9610 



! 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Sphere2(Dt( ♦ ) ,J ,Tc ,Patn,Rh,Tgttn,Dia ,Efreq ,Nuspkr .Diaspkr.Gpre ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Er WITH 
REM TARGET NO. 2 EMPLOYED. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++ LOADING DATA +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

f 

! 

J=39 

RESTORE S2 
FOR 1=0 TO J 
READ Dt( I ) 

NEXT I 

READ Tc ,Patn,Rh,Tgttn,Dia ,Efreq ,Nuspkr,Diaspkr,Gpre 
2: DATA 4.72,4.66.4.16,3.23,2.01 , . 560 . 940 .-2 . 35 . -3 . 52 ,~4 . 36 

DATA -4. 77,-4. 70, ~4. 18,-3.25,-! .99, -.560,. 960, 2. 34. 3. 51 ,4.34 
DATA 4.74 ,4.66 ,4. 13 ,3. 1 9. 1 .95,. 500.-1 .00,-2.40,-3.57,-4.38 
DATA -4. 78, -4. 68, -4. 15,-5. 20,-1 .92, -.480,1 .02,2.39,3.54,4.34 
DATA 21 .0, 1009. 4 ,49.8 ,30.55 . . 1013 ,5000 . 1 9 , .0762 ,11094 

i 

I 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

I++++ MODULE CONCLUSION ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

1 

I + ++++++++ +++ + + + + + + + + + + +++ + + ++++++ + +++ + + +++ + +++ + + ++++++++++ +++++ +++ ++++-♦- 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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9620 ! 

9630 ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

9640 SUB Sphere3(Dt(*),J ,Tc,Patrn,Rh,Tgttm,Dia,EfreQ,Nu5pkr,Diaspkr,Gpre) 

9650 REM 4-++++++++++++++++++++++++^+++++ +++++++++++ ++++++++++++++++++++++++ 

9660 REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Er WITH 

9670 REM TARGET NO. 3 EMPLOYED. 

96 B0 REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9690 ! 

9700 !+-»- + + LOADING DATA + + + + + + + + + +++ + +++++ + + ++++ + + + + + + +++ + +++ ++++++ + + +++++ +++ 

9710 I 

9720 ! 

9730 J=39 

9740 RESTORE S3 
9750 FOR 1=0 TO J 
9760 READ DU I ) 

9770 NEXT I 

9780 READ Tc ,Patm ,Rh ,Tgt tn ,Dia ,Efreq .Nuspkr .Diaspkr ,Gpre 
9790 S3: DATA 3 . 6 1 . 3 . 40 ,2 . 86 ,2 . 04 . 1 . 0 1 ) 20 1 . 25 ,-2 . 25 ,-3 . 04 .-3 . 54 

9800 DATA -3.68 ,-3.46 ,-2.9) ,-2.07 ,-l .04 ,. 110 .1 .25 ,2.24 ,3.00 ,5.48 

9810 DATA 3.63 ,3.39,2.84,2.02 ,.970,-. 180,-1 .50,-2.51 ,-3.08,-3.56 

9820 DATA -3 . 6 9 , -3 . 46 , -2 . 89 , -2 . 04 . 990 , . 1 60 , 1 . 29 , 2 . 26 , 3 . 05 ,5 . 50 

9830 DATA 2 1 . 3 , 1 009 . 6 . 49 . 8 ,30 . 59 . .07634 ,5000 , 1 9 , .0762 ,11094 

9840 ! 

9850 ! 

9860 l+ + 4-4-++4-++ + 4- + + ++ + + + + + + + + + +++ + + +'^ + + + -f + 4’ + + + -‘- + + + + + + + + + + + +-l- + + + + + + + + + + + + + + + + + 

9870 1 

9880 I++++ MODULE CONCLUSION +-f*+4-+++ + + +++4'+ + +++++ + + +++ + + + + +++ + ++++++ + + +++ +++ + 

9890 ! 

9900 SUBEND 

9910 ! 

9920 !+ + + + + + + + + + + + + + + 4-4- + + + + + + + 4- + + -l- + + + + + + + + -f- + + + + + + 4--l> + -f + +'+ + + + 4--l- + + + 4- + + + + + + + + + + + + 

9930 ! 

9940 I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
9950 ! 
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9960 
9970 
9980 
9990 
1 0000 
10010 
10020 
10030 
10040 
10050 
10060 
1 0070 
10080 
10090 
10100 
10110 
10120 
10130 
10140 
10150 
10160 
10170 
10180 
10190 
10200 
10210 
10220 
10230 
10240 
10250 
10260 
10270 
10280 
10290 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB Sphere4(Dt< ♦ ) ,J ,Tc ,Patn ,Rh,Tgttn,D:a ,Efreq ,Nuspkr,Diaspkr ,Gpre ) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CONTAINS THE DATA USED FOR DETERMINATION OF Er WITH 
REM TARGET NO. 4 EMPLOYED. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 



!+ +++ LOADING DATA +++ + -I 



h++++++++++H 



^++++++++ 



I 



J=39 

RESTORE S4 
FOR 1=0 TO J 
READ Dt(I ) 

NEXT I 

READ Tc,Patn,Rh,Tgttm,Dia ,Efreq,Nuspkr,Diaspkr,Gpre 
S4: DATA 2 . 28 .2 . 1 0 , 1 . 68 . 1 . 1 0 . . 420 . 320 1 . 03 1 . 64 .-2 . 07 .-2 . 3 1 

DATA -2.34 , -2. 12 .-1 .71 .-1 . 12 ,-.440 . .320 J .02 , 1 .62 ,2.06 ,2.28 
DATA 2.29,2. 10 , 1 .68 , 1 .09 ,.400 ,-.340 ,-1..04 ,-l .66 ,-2. 10,-2.32 
DATA -2. 35,-2. 12 ,-l .70 .-1 . 1 1 ,-.400 . .340 . 1 .04 j .B4 ,2.08 ,2.28 
DATA 21 .4 ,1009.9,4 9.6 ,30.64 ,.06240, 5000, 19,. 0762 ,1 1094 



I 



++++++++++++H 



h + + + H 



++++ MODULE CONCLUSION -^+ + -^h 
SUBEND 

++++++++++++++++ +++++++++H 



h + + + + + + -^ + + H 



h++++++++++ 



+ + ^ ++ + + + + + + + + + 



I 



++++++++++++++++++++++++++++++++++ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$s$$$$$$$$$$ $$$$$$$$$$ 
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10300 
10310 
10320 
10330 
10340 
10350 
10360 
10370 
1 0380 
10390 
1 0400 
10410 
10420 
10430 
1 0440 
10450 
10460 
10470 
10480 
1 0490 
10500 
10510 
10520 
10530 
10540 
10550 
10560 
10570 
10580 
10590 
10600 
10610 



!$$S$$$S$$$$S$$SS$$SS$$S$$$$S$$$S$S$$$$$S$$SSS$$$$$$S$$S$S$S$$$$$$$$SS$$ 
sue Trans30( J ,Dt ( * ) ) 

REM + + + + + ++ + + +++ + + + + + + + +++ + ++-I- + + +++++ + + + + + + + + + + + + + + + + +++ + + + + +++ + + + + + + + 

REM THIS MODULE CONTAINS DATA USED IN DETERMINATION OF Er WITH 

REM 3.0U APPLIED. 

REM +++++++ + +++++++++++++++++++++++++++++++-I-+-I- + +++ +++++++++++++++++++ + 

j 

I 

j 

J=39 

RESTORE T30 
FOR 1=0 TO J 
READ Dt( I ) 

NEXT I 

T30: DATA . 6 90 . . 680 . . 590 , . 440 , . 520 , . 1 60 , - . 020 , - . 230 . - . 4 1 0 . 550 

DATA -.61 0 ,-.610 ,-.520 ,*-.360 ,-.240 , - . 080 , . 1 1 0 , . 320 , . 490 . .630 
DATA .690 . .680 , . 590 . . 430 , . 320 , . 1 50 , - . 040 ,- . 240 . 420 . 560 

DATA -.620 ,-.600 ,-.510 ,-.380 ,-.240 ,-.070 .. 120 , .330 . .500 ,.640 



+++++++++++++++H 



^ + + + + -l- + + - 



h + + + 



++-I-+ MODULE CONCLUSION ++++++++ ++++++++++++++++ + H 

j 

SUBEND 

I 

++++++++++++++++++++++++++++r++++++++++++++++++++++++++ +++++++ +++++++++ 

s $ s$$$ S $ $$SS $ S $$$$ $$S SS $S S S SS SSSS$$$SSSSSSSS$SS$S$SSS $$$£$£$$$$ $$ss$s$s 
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10G20 
10G30 
10G40 
10G50 
10GG0 
10G70 
10G80 
I0G90 
10700 
10710 
10720 
10750 
10740 
107B0 
107G0 
10770 
10780 
10790 
1 0800 
10810 
10820 
10830 
10840 
10850 
108G0 
10870 
10880 
10890 
10900 
10910 
10920 
10950 



! 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Trans35( J , Dt ( ♦ ) ) 

REM +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ + 

REM THIS MODULE CONTAINS DATA USED IN DETERMINATION OF Er WITH 
REM 3.5U APPLIED. 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++H 



!++++ LOADING DATA +++++++++++++H 



I- + + + + + + + + + ++ + + + + + + H 



J-39 

RESTORE T35 
FOR 1=0 TO J 
READ Dt( I ) 

NEXT I 

T35: DATA . 830 , . 770 , . G40 . . 490 , . 290 , . 080 . 1 80 . 400 . - . G 1 0 . 730 

DATA -.770 ,-.720 ,-.580 ,-.430 ,- . 240 , - . 020 , . 240 , .4G0, .GG0 , .780 
DATA .850 , .770 . . G30 , . 470 . . 290 , . 0G0 . 1 90 ,-.420 ,-.G00 ,-.720 
DATA - . 770 . - . 7 1 0 . - . 570 . - . 420 , - . 230 , - . 0 1 0 , . 240 . . 480 , . G70 , . 790 

! 



h + + + H 



I- + + + + + + + + 



h++++++++ 



I 



++++++++++++++H 

++++ MODULE CONCLUSION +++++H 
SUBEND 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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10940 
10950 
10960 
10970 
10980 
10990 
1 1 000 
11010 
1 1020 
1 1030 
1 1040 
1 1050 
1 1 060 
1 1070 
1 1080 
1 1090 
11100 
11110 



I 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Noise 30 (J,Jt ,Dt( * ),Tdt( ♦ ),Tcn,Rhn,Patnn,Efreqn) 

REM ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM THIS MODULE CONTAINS DATA USED IN DETERMINATION OF Er WITH 
REM 3.0U APPLIED. 

! 

!++++ LOADING DATA + 



I- + + + + H 
‘ I 
I 



J = 59 

RESTORE N30 
FOR 1=0 TO J 
READ Dt( I ) 
NEXT I 
Jt=39 

RESTORE Tn50 



11120 


FOR 


1=0 TO 


Jt 






11130 


READ Tdt( I ) 






11140 


NEXT 


I 








11150 


READ 


Ten ,Rhn ,Patnn ,Efreqn 






11160 


N30: 


DATA 


0,0, -.02 ,-.03,-.04,- 


.06,- 


.08 ,-. 1 ,-. 1 1 ,-. 12 


11170 




DATA 


12 1 1 1 ,-.08 


.07,- 


.05 ,-.04 ,-.01 ,0,0 


11180 




DATA 


0,0,-.02,-.03,-.05,- 


.07,- 


.08,-. 1 ,-. 11 ,-. 12 


11190 




' DATA 


12 12 1 1 ,-.08 . 


-.07 , 


-.04 ,-.03 ,-.02 ,0 ,0 


1 1 200 




DATA 


.01 ,0 ,0 ,-.03 ,-.04 


07 ,-. 


09 1 1 ,- . 12 ,- . 1 2 


11210 




DATA 


-. 12 ,-. 12 ,-. 1 1 ,-.08 . 


-.07 , 


- . 04 .05 ,0 ,0 , . 01 


1 1220 


T n50 


: DATA 


.69, .68, .59, .44, .32, 


. IB ,- 


.02 ,-.23 ,-.41 ,-.55 


1 1230 




DATA 


-.61 ,-.61 ,-.52 ,-'.36 , 


-.24, 


-.08 , . 11 , .32 . .49 , .! 


1 1 240 




DATA 


.69 . .68 , .59 , .45 , .32 , 


.15,- 


.04 ,-.24 .-.42 ,-.56 


1 1250 




DATA 


-.62 ,-.B ,-.51 ,-.38,- 


.24,- 


.07, .12, .33, .5, .64 


1 1260 




DATA 


20.8,51 .7,1016.7,5000 





1 1270 
1 1280 
1 1290 
1 1300 
11310 
1 1320 
1 1 330 
1 1 340 
1 1350 
1 1360 
1 1 370 
1 1380 



K++++++++++++++++H 






I- + + 



I- + + + J- + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + 
! 



I +++++++++++H 
I 

!++++ MODULE CONCLUSION + + + ++-^-i 
SUBEND 

I+++++++++++++++++++++++++++H 
SSSSSSSSSSSS%SSSSSSSSSSSSSSSSSS$$SSSSSSSSSSSSSSSSSSSSSSSSS$$SSSSSSSSSSS 



I 



h ++++++++++++++++++++ ++++++ 
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1 1390 
1 1 400 
11410 
1 1 420 
1 1430 
1 1440 
1 1 450 
1 1 4B0 
1 1 470 
1 1480 
1 1490 
1 1500 
1 1510 
1 1520 
1 1530 
1 1540 
1 1550 
1 1560 
1 1570 
1 1580 
1 1590 
1 1600 
11610 
1 1620 
1 1650 
1 1640 
1 1650 
1 1660 
1 1670 
1 1680 
1 1690 
1 1700 
1 1710 
1 1720 
1 1730 
1 1740 
1 1750 
1 1760 
1 1 770 
1 1780 
1 1790 
1 1 800 
11810 



I 



I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB Noise 35 ( J ,Jt ,Dt( * ) ,Tdt( ♦ ) ,Tcn,Rhn,Patnn,Efreqn ) 

REN + + ++++ + + +++++++ + +++-I- + + +++++++++ + + +++++++++++++ +++++ + + +++ ++++++++++ 

REN THIS NODULE CONTAINS DATA USED IN DETERMINATION OF Er WITH 
REM 3.50 APPLIED. 

REM ++++++-}■ + + +++ + + + ++++++++++++++++++ ++++++++++++-}■+++ ++++ + + + +++ +++ + + + + 

! 

!+ + + + LOADING DATA ++ + +++++ ++++++++++ + + ++++++ ++++++ + +++++ + -I 

I 



J = 39 

RESTORE N55 
FOR 1=0 TO J 
READ Dtd ) 

NEXT I 
Jt = 39 

RESTORE Tn35 
FOR 1=0 TO Jt 
READ Tdt( I ) 

NEXT I 

READ Ten ,Rhn ,Pat nn ,Ef reqn 

N55: DATA - . 06 . 08 . 08 . 08 . 07 . 06 . 04 . 02 .0 . . 02 

DATA . 05 . . 06 , . 06 . . 06 . . 04 . . 03 , . 02 .0 , - . 03 , - . 05 
DATA -.07,-.08 ,-.09.-.08 ,-.08,-.06 ,-.04.-.02 ,.01 ,.03 
DATA .04,. 07.. 06,. 07,. 05,. 04,. 01 ,-.01 ,-.04,-. 06 
Tn35: DATA . 74 , . 83 , . 82 , . 72 , . 53 , . 38 , . 1 8 ,~ . 06 ,- . 31 ,- . 53 
DATA - . 68 , - . 77 , - . 75 ,-.64, -.46, -.3,-.! ,.14, .4, .6 
DATA .75 , .83 , .81 ,.7, .5 . .36 ,. 16 ,-.08 ,-.33 ,-.54 
DATA -.69 ,-.76 .74 ,-.63 ,-.46 .-.28 ,-.08 ,. 16 , .41 , .62 
DATA 21 .3,51 . 3 , 1 0 1 6 . 6 .5000 



! 



I +++++++++++++++++++++++++++++++++++++++++++ 
I 



l- + + + + + ^ 



- + + 



!+ + + + MODULE CONCLUSION + + + +++ + -I-+++++ + 

I 




+++++++++++++++++ 



SUBEND 



j 

I +++++++++T+++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++ 
! 



! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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APPFNmX F 



ACOUSTIC ARRAY CALIBRATION DATA 

This appendix contains data used to determine the calibration 
product E^. Data units are contained in square brackets. 

For listings of the same type of data the units and uncertainty are 
explicitly given for the first datum. Remaining data in the listing have the 
same units and uncertainty as the first datum. 

1 . TWO WAY PROPAGATION PATH DATA 

Pre-amplifier gain and array impedance 2 were previously 

determined by Moxcey. G^was determined to be 1 1,094 Z was determined 
to be 15.8[0]. [Ref. 8] 

The remaining data are presented in the following tables. 

'^supply voltage supplied by the HP3314A function generator to 

the pre-amplifier (see Figure 8). 



MO 



Table 5a: Two Way Propagation Path Data for Sphere One, 
Collected 30 November 1987; Run *' 2 \ Data Re-measured 
Because First Run (with 3.5 [V] Supply) Was Clipped. 



Diameter = 2a = .2546 + 0.0006 [m] 

T^ = 20.2 ±0.1 [° c] 

Rh = 52.8 ± 0.1 [%] 

P = 1016.8 ± 0.2 [mb] 

No. pulses = 20 

%piy = 3.0 [V] 

t^ = 30.3 1 [mS] - 400 [pS] = 29.9 1 ± 0. 1 [mS] 



t 


Vt 


2.25 ± O.OUmS] 


690 ± 


2.26 


680 


2.27 


590 


2.28 


440 


2.29 


320 


2.30 


160 


2.31 


- 20 


2.32 


-230 


2.33 


-410 


2.34 


-550 


2.35 


-610 


2.36 


-610 


2.37 


-520 


2.38 


-360 


2.39 


-240 


2.40 


- 80 


2.41 


1 10 


2.42 


320 


2.43 


490 


2.44 


630 



t 

5 [mV] 2.45 

2.46 

2.47 

2.48 

2.49 

2.50 

2.51 

2.52 

2.53 

2.54 

2.55 

2.56 

2.57 

2.58 

2.59 

2.60 

2.61 

2.62 

2.63 

2.64 



Vt 

690 

680 

590 

430 

320 

150 

- 40 
-240 
-420 
-560 
-620 
-600 
-510 
-380 
-240 

- 70 
120 
330 
500 
640 
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Table 5b; Two Way Propagation Path Data for Sphere One Continued, 
Collected 30 November 1987; Run *2; Data Re-measured 
Because First Run (with 3.5 [V] Supply) Was Clipped. 



t 

33.68 lO.OUmS] 

33.69 

33.70 

33.71 

33.72 

33.73 

33.74 

33.75 

33.76 

33.77 

33.78 

33.79 

33.80 

33.81 

33.82 

33.83 

33.84 

33.85 

33.86 

33.87 



9.54 ± 0.03 [V] 

9.40 

8.32 

6.42 

3.88 
0.960 
-2.07 
-4.90 
-7.24 
-8.87 
-9.64 
-9.45 
-8.32 
-6.41 
-3.83 
-0.890 

2.14 

4.94 

7.26 

8.88 



t 


V 

r 


33.88 


9.60 


33.89 


9.38 


33.90 


8.22 


33.91 


6.26 


33.92 


3.67 


33.93 


0.730 


33.94 


-2.30 


33.95 


-5.1 1 


33.96 


-7.40 


33.97 


-8.96 


33.98 


-9.64 


33.99 


-9.39 


34.00 


-8.20 


34.01 


-6.21 


34.02 


-3.61 


34.03 


-0.650 


34.04 


2.36 


34.05 


5.12 


34.06 


7.38 


34.07 


8.91 
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Table 6a: Two Way Propagation Path Data for Sphere Two, 
Collected 28 November 1987. 



Diameter = 2a = 1 0. 1 3 ± 0.0 1 (cm] 

T^ = 21.0 ± 0.1 [°c] 

Rh = 49.8 ± 0.1 [%] 

P = 1009.4 + 0.2 [mb] 

No. pulses = 20 

%p.y = ^-5 IV] 

t^ = 30.95 [mS] - 400 [pS] = 30.55 + 0. 1 [mS] 



t 

2.26 + 0.0 KmS] 

2.27 

2.28 

2.29 

2.30 

2.31 

2.32 

2.33 

2.34 

2.35 

2.36 

2.37 

2.38 

2.39 

2.40 

2.41 

2.42 

2.43 

2.44 

2.45 



830 + 5 (mV] 
770 
640 
490 
290 
80 
-180 
-400 
-610 
-730 
-770 
-720 
-580 
-430 
-240 
- 20 
240 
460 
660 
780 



2.46 


830 


2.47 


770 


2.48 


630 


2.49 


470 


2.50 


290 


2.51 


60 


2.52 


-190 


2.53 


-420 


2.54 


-600 


2.55 


-720 


2.56 


-770 


2.57 


-710 


2.58 


-570 


2.59 


-420 


2.60 


-230 


2.61 


- 10 


2.62 


240 


2.63 


480 


2.64 


670 


2.65 


790 
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Table 6b: Two Way Propagation Path Data for Sphere Two Continued, 
Collected 28 November 1987. 



t 


Vr 


t 


V 

r 


34.13 lO.OHmS] 


4.72 + 0.01 [V] 


3433 


474 


34.14 


4.66 


3434 


466 


34.15 


4.16 


3435 


413 


34.16 


3.23 


3436 


3.19 


34.17 


2.01 


3437 


1.95 


34.18 


0.560 


3438 


0.500 


34.19 


-0.940 


3439 


-1.00 


34.20 


-2.35 


3440 


-2.40 


34.21 


-3.52 


3441 


-3.57 


34.22 


-4.36 


3442 


-438 


34.23 


-4.77 


34.43 


-478 


34.24 


-4.70 


34.44 


-468 


34.25 


-418 


3445 


-415 


34.26 


-3.25 


34.46 


-3.20 


34.27 


-1.99 


3447 


-1.92 


34.28 


-0.560 


34.48 


-0.480 


34.29 


0.960 


3449 


1.02 


34.30 


2.34 


3450 


2.39 


34.31 


3.51 


3451 


3.54 


34.32 


434 


34.52 


434 
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Table 7: Two Way Propagation Path Data for Sphere Three, 
Collected 28 November 1987. 



Diameter = 2a = 7.634 + 0.003 [cm] 

T^ = 21.3 i O.I P c] 

Rh = 49.8 ± 0. 1 [%] 

P= 1009.6 + 0.2 [mb] 

No. pulses = 20 
%p,y ' 35 [VI 

t,. - 30.99 ImSI - 400 1^51 = 30.59 1 0. 1 [mSl 
data is the same as sphere two's (see Table 6a). 



t 


Vp 


t 


Vp 


33.98 ±0.01[mS] 


3.61 ± 0.03 [V] 


34.18 


3.63 


33.99 


3.40 


34.19 


3.39 


34.00 


2.86 


34.20 


2.84 


34.01 


2.04 


34.21 


2.02 


34.02 


1.01 


34.22 


0.970 


34.03 


-0.120 


34.23 


-0.180 


34.04 


-1.25 


34.24 


-1.30 


34.05 


-2.25 


34.25 


-2.31 


34.06 


-3.04 


34.26 


-3.08 


34.07 


-3.54 


34.27 


-3.56 


34.08 


-3.68 


34.28 


-3.69 


34.09 


-3.46 


34.29 


-3.46 


34.10 


-2.91 


34.30 


-2.89 


34. 11 


-2.07 


34.31 


-2.04 


34.12 


-1.04 


34.32 


-0.990 


34.13 


0.1 10 


34.33 


0.160 


34.14 


1.25 


34.34 


1.29 


34.15 


2.24 


34.35 


2.26 


34.16 


3.00 


34.36 


3.05 


34.17 


3.48 


34.37 


3.50 
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Table 8; Two Way Propagation Path Data for Sphere Four, 
Collected 28 November 1987. 



Diameter = 2a = 6.240 + 0.008 [cm] 

T^ = 21.4 ± 0.1 [° c] 

Rh = 49.6 ± 0. 1 [%] 

P= 1 009.9 ± 0.2 [mb] 

No. pulses = 20 
^supply = 3.5 [VI 

t, = 3 1 .04 ImS] - 400 luS) = 30.64 ♦ 0. 1 ImSl 
data is the same as sphere two's (see Table 6a). 



t 


V 

r 


t 


Vr 


33.83 ±0.0 UmS] 


2.28 ± 0.03 [V] 


34.03 


2.29 


33.84 


2.10 


34.04 


2.10 


33.85 


1.68 


34.05 


1.68 


33.86 


1.10 


34.06 


1.09 


33.87 


0.420 


34.07 


0.400 


33.88 


-0.320 


34.08 


-0.340 


33.89 


-1.03 


34.09 


-1.04 


33.90 


-1.64 


34.10 


-1.66 


33.91 


-2.07 


34. 11 


-2.10 


33.92 


-2.31 


34.12 


-2.32 


33.93 


-2.34 


34.13 


-2.35 


33.94 


-2.12 


34.14 


-2.12 


33.95 


-1.71 


34.15 


-1.70 


33.96 


-1.12 


34.16 


-1.11 


33.97 


-0.440 


34.17 


-0.400 


33.98 


0.320 


34.18 


0.340 


33.99 


1.02 


34 19 


1.04 


34.00 


1.62 


34.20 


1.64 


34.01 


2.06 


34.21 


2.08 


34.02 


2.28 


34.22 


2.28 
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Table 9: Noise Data for 3.0 [V] Supply Voltage, 
Collected 30 November 1987. 



= 20.8 1 0. 1 [° c] Rh = 51.7 i 0.1 (%] P = 1016.7 ± 0.2 [mb] 



No. pulses = 20 



data is the same as sphere one's (see Table 5a). 




t 


Vn 


t 


Vn 


33.73 iO.OUmS] 


00 i 5 [mV] 


34.03 


-120 


33.74 


00 


34.04 


-120 


33.75 


- 20 


34.05 


-1 10 


33.76 


- 30 


34.06 


- 80 


33.77 


- 40 


34.07 


- 70 


33.78 


- 60 


34.08 


- 40 


33.79 


- 80 


34.09 


- 30 


33.80 


-100 


34.10 


- 20 


33.81 


-1 10 


34.1 1 


00 


33.82 


-120 


34.12 


00 


33.83 


-120 


34.13 


10 


33.84 


-1 10 


34.14 


00 


33.85 


-100 


34.15 


00 


33.86 


- 80 


34.16 


- 30 


33.87 


- 70 


34.17 


- 40 


33.88 


- 50 


34.18 


- 70 


33.89 


- 40 


34.19 


- 90 


33.90 


- 10 


34.20 


-1 10 


33.91 


00 


34.21 


-120 


33.92 


00 


34.22 


-120 


33.93 


00 


34.23 


-120 


33.94 


00 


34.24 


-120 


33.95 


- 20 


34.25 


-1 10 


33.96 


- 30 


34.26 


- 80 


33.97 


- 50 


34.27 


- 70 


33.98 


- 70 


34.28 


- 40 


33.99 


- 80 


34.29 


- 30 


34.00 


-100 


34.30 


00 


34 01 


-1 10 


34.31 


00 


34.02 


-120 


34.32 


10 
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Table 10a; Noise Data for 3.5 [V] Supply Voltage, 
Collected 30 November 1987. 



T^ = 21.3 ±0.1 [°c] 

Rh = 51.3 i 0.1 [%] 

P = 1016.6 1 0.2 [mb] 
No. pulses = 20 



t 


Vt 


t 


Vt 


2.24 ±0.0 l[mS] 


740 ± 5 [mV] 


2.44 


750 


2.25 


830 


2.45 


830 


2.26 


820 


2.46 


810 


2.27 


720 


2.47 


700 


2.28 


530 


2.48 


500 


2.29 


380 


2.49 


360 


2.30 


180 


2.50 


160 


2.31 


- 60 


2.51 


- 80 


2.32 


-310 


2.52 


-330 


2.33 


-530 


2.53 


-540 


2.34 


-680 


2.54 


-690 


2.35 


-770 


2.55 


-760 


2.36 


-750 


2.56 


-740 


2.37 


-640 


2.57 


-630 


2.38 


-460 


2.58 


-460 


2.39 


-300 


2.59 


-280 


2.40 


-100 


2.60 


- 80 


2.41 


140 


2.61 


160 


2.42 


400 


2.62 


410 


2.43 


600 


2.63 


620 



148 



Table 10b; Noise Data for 3.5 [V] Supply Voltage Continued; 
Collected 30 November 1987. 



t 


Vn 


t 




33.80 ±0.01[mS] 


-60 + 5 [mV] 


34.00 


- 70 


33.81 


- 80 


34.01 


- 80 


33.82 


- 80 


34.02 


- 90 


33.83 


- 80 


34.03 


- 80 


33.84 


- 70 


34.04 


- 80 


33.85 


- 60 


34.05 


- 60 


33.86 


- 40 


34.06 


- 40 


33.87 


- 20 


34.07 


- 20 


33.88 


00 


34.08 


10 


33.89 


20 


34.09 


30 


33.90 


50 


34.10 


40 


33.91 


60 


34.1 1 


70 


33.92 


60 


34.12 


60 


33.93 


60 


34.13 


70 


33.94 


40 


34.14 


50 


33.95 


30 


34.15 


40 


33.96 


20 


34.16 


10 


33.97 


00 


34.17 


- 10 


33.98 


- 30 


34.18 


- 40 


33.99 


- 50 


34.19 


- 60 
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2. ONE WAY PROPAGATION PATH DATA 



One way propagation path data are presented in the following tables. 
T^ and Rh are estimated. The estimates were based on hygrothermographic 

comparisons and 28 and 30 November data. 



Table 1 1; One Way Propagation Path Data for 3.0 [V], 
5000 [Hz] Continuous Wave Supply from the HP33 1 4A 
Function Generator; Collected 7 December 1987. 



(Vrn,5\ = 0.246 1 0.001 [V] 

T^~ 21.0 [° c] 

Rh ~ 50.0 [%] 

P= 1021.5 + 0.2 [mb] 

R = 5.30 [m] (calculated from two way propagation path data for 
spheres two, three and four) 

I (referenced to 10”’^ [W/m^]) 

r 

94.7 + 0. 1 [dB] 

94.5 

94.7 

95.0 

94.8 

95.1 

94.2 
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Table 12; One Way Propagation Path Data for 3.5 [V], 
5000 [Hz] Continuous Wave Supply from the HP33 1 4A 
Function Generator; Collected 7 December 1987. 



= 0.523 ±0.001 [V] 

T^~21.0[° c] 

Rh « 50.0 [%] 

P = 1021.5 ± 0.2 [mb] 

R = 5.30 [m] (calculated from two way propagation path data for 
spheres two, three and four) 



1 (referenced to 10 [W/m^]) 

r 

97.4 ± 0. 1 [dB] 

96.4 

95.1 

97.0 

96.2 

97.1 

96.4 

95.4 

96.3 

97.4 
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